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Abstract 



We review various approaches to the calculation of QCD condensates and of the nucleon 
characteristics in nuclear matter. We show the importance of their self-consistent treatment. 
The first steps in such treatment appeared to be very instructive. It is shown that the alleged 
pion condensation anyway can not take place earlier than the restoration of the chiral symmetry. 
We demonstrate how the finite density QCD sum rules for nucleons work and advocate their 
possible role in providing an additional bridge between the condensate and hadron physics. 
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1 Introduction 



The nuclear matter, i.e. the infinite system of interacting nucleons was introduced in order to simplify 
the problem of investigation of finite nuclei. By introducing the nuclear matter the problems of NN 
interaction in medium with non-zero baryon density and those of individual features of specific 
nuclei were separated. However, the problem of the nuclear matter is far from being solved. As we 
understand now, it cannot be solved in consistent way, being based on conception of NN interactions 
only. This is because the short distances, where we cannot help considering nucleons as composite 
particles, are very important. 

There is limited data on the in-medium values of nucleon parameters. These are the quenching 
of the nucleon mass m and of the axial coupling constant qa at the saturation value po with respect 
to their vacuum values. The very fact of existing of the saturation point po is also the "experimental 
data", which is the characteristics of the matter as a whole. The nowadays models succeeded in 
reproducing the phenomena although the quantitative results differ very often. 

On the other hand, the knowledge about the evolution of hadron parameters is important for 
understanding the evolution of the medium as a whole while the density p of distribution of the 
baryon charge number increases. (When p is small enough, it is just the density of the distribution 
of nucleons). There can be numerous phase transitions. At certain value of density p = p a the 
Fermi momenta of the nucleons will be so large, that it will be energetically favourable to increase 
p by adding heavier baryons instead of new nucleons. The nuclear or, more generally, hadronic 
matter may accumulate excitations with the pion quantum numbers, known as pion (or even kaon) 
condensations. Also the matter can transform to the mixture of hadrons and quark-gluon phase or 
totally to the quark-gluon plasma, converting thus to baryon matter. The last but not the least is the 
chiral phase transition. The chiral invariance is assumed to be one of the fundamental symmetries 
of the strong interactions. 

The chiral invariance means that the Lagrangian as well as the characteristics of the system 
are not altered by the transformation ip — > ipe 10115 of the fermion fields ip. The model, suggested by 
Nambu and Jona-Lasinio (NJL) provides a well-known example. The model describes the massless 
fermions with the four-particle interactions. In the simplest version of NJL model the Lagrangian is 



If the coupling constant G is large enough, the mathematical (empty) vacuum is not the ground 
state of the system. Due to the strong four-fermion interaction in the Dirac sea the minimum of 
the energy of the system is reached at a nonzero value of the fermion density. This is the physical 
vacuum corresponding to the expectation value (Ol^^O) 7^ 0. 

This phenomenon is called "spontaneous chiral symmetry breaking". In the physical vacuum the 
fermion obtains the mass 



caused by the interaction with the condensate. On the other hand, the expectation value (0\tpip\0) 
is expressed through the integral over the Dirac sea of the fermions. Of course, we have to introduce 




(1) 



m = -2G(0\ifnl>\0) 



(2) 
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a cutoff A to prevent the ultraviolet divergence caused by the four-fermion interaction 



m f 



A 



(01-0-010) 




dp 



(p 2 + m 2 ) 1 / 2 



(3) 



Thus Eqs.(2) and (3) compose self-consistent set of equations which determine the values of the 
condensate and of fermion mass m in the physical vacuum. 

Originally the NJL model was suggested for the description of the nucleons. Nowadays it is used 
for the quarks. The quark in the mathematical vacuum, having either vanishing or very small mass 
is called the "current" quark. The quark which obtained the mass, following Eq.(2) is called the 
"constituent" quark. In the nonrelativistic quark model the nucleon consists of three constituent 
quarks only. 

Return to the nuclear matter. To understand, which of the hadron parameters are important, note 
that we believe nowadays most of the strong interaction phenomena at low and intermediate energies 
to be described by using effective low-energy pion-nucleon or pion-constituent quark Lagrangians. 
The ttN coupling constant is: 

= 9± (4) 

with f n being the pion decay constant. This is the well-known Goldberger- Treiman (GT) relation 
@. It means, that the neutron beta decay can be viewed as successive strong decay of neutron to 
tt~P system and the decay of the pion. Thus, except the nuclear mass m*(p), the most important 
parameters will be the in- medium values g\(p), fZ(p) an d m l(p)- 

On the other hand, the baryonic matter as a whole is characterised by the values of the con- 
densates, i.e. by the expectation values of quark and gluon operators. Even at p = some of the 
condensates do not vanish, due to the complicated structure of QCD vacuum. The nonzero value of 
the scalar quark condensate (Ol^^O) reflects the violation of the chiral symmetry. In the exact chiral 
limit, when (O^^O) = (and the current quark masses vanish also), the nucleon mass vanishes too. 
Thus, it is reasonable to think about the effective nucleon mass m*(p) and about the other parame- 
ters as the functions of the condensates. Of course, the values of the condensates change in medium. 
Also, some condensates which vanish in vacuum may have the nonzero value at finite density. 

At the same time, while calculating the expectation value of the quark operator ipip in medium, one 
finds that the contribution of the pion cloud depends on the in-medium values of hadron parameters. 
Hence, the parameters depend on condensates and vice versa. Thus we came to the idea of self- 
consistent calculation of hadron parameters and of the values of condensates in medium. The idea 
of self-consistency is, of course, not a new one. We have seen just now, how NJL model provides 
an example. We shall try to apply the self-consistent approach to the analysis of more complicated 
systems. 

The paper is organized as follows. In Sec. 2 we review the present knowledge on the in-medium 
condensates. In Sec. 3 we present the ideas and results of various approaches to calculation of the 
hadron parameters in medium. We review briefly the possible saturation mechanisms provided by 
these models. In Sec. 4 we consider the first steps to self-consistent calculation of scalar condensate 
and hadron parameters. The experience appeared to be very instructive. Say, the analysis led to the 
conclusion, that in any case the chiral phase transition takes place at the smaller values of density 
than the pion condensation. Hence, the Goldstone pions never condense. However, analysis of the 
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behaviour of the solutions of the corresponding dispersion equation at larger densities appears to be 
useful. 

Suggesting QCD sum rules at finite density as a tool for a future complete self-consistent investi- 
gation, we show first how the method works. This is done in Sec. 5. In Sec. 6 we present more detailed 
self-consistent scenario. 

We present the results for symmetric matter, with equal densities of protons and neutrons. 

Everywhere through the paper we denote quark field of the flavour % and colour a as if>£. We 
shall omit the colour indices in most of the cases, having in mind averaging over the colours for 
colourless objects. As usually, <jj, Tj and 7 M are spin and isospin Pauli matrices and 4 Dirac matrices 
correspondingly. For any four-vector we denote A fl / y fl = A^j^ = A. The system of units with 
h — c — 1 is used. 



2 Condensates in nuclear matter 
2.1 Lowest order condensates in vacuum 

The quark scalar operator ipif) is the only operator, containing minimal number of the field operators 
x/j, for which the expectation value, in vacuum has a nonzero value. One can find in the textbooks 
a remarkable relation, based on partial conservation of axial current (PCAC) and on the soft-pion 
theorems 

mlJl = -i(0| [^(0)^(0)^(0)]] |0> (5) 

with m^b, f n standing for the mass and the decay constant of pion, H being the density of the 
Hamiltonian of the system, while F b 5 are the charge operators, corresponding to the axial currents, 
b is the isospin index. 

Presenting (effective) Hamiltonian 

H = H + H b (6) 

with Ho(Hb) conserving (explicitly breaking) the chiral symmetry, one finds that only Hb piece 
contributes to Eq.(4). In pure QCD 

H b = H® CD = m u uu + m d dd , (7) 

with m u ^ standing for the current quark masses. This leads to well known Gell-Mann-Oakes-Renner 
relation (GMOR) |] 

2f 2 m 2 

(0\uu + dd\0) = . (8) 

m u + m d 

Of course, assuming SU(2) symmetry, which is true with the high accuracy, one finds (0|mt|0) = 
(0|dd|0). Numerical value (0|uu|0) = (-240 MeV) 3 can be obtained from Eq.@. 

The quark masses can be obtained from the hadron spectroscopy relations and from QCD sum 
rules — see the review of Gasser and Leutwyler M. Thus the value of the quark condensate was 
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calculated by using Eq.(^). The data on the lowest order gluon condensate (a is the colour index, a s 
is the QCD coupling constant) 

(0| ^ G^GT |0) « (0.33 GeV) 4 (9) 

was extracted by Vainshtein et al. || from the analysis of leptonic decays of p and <p mesons and 
from QCD sum- rules analysis of charmonium spectrum JJ. 



2.2 Gas approximation 

In this approximation the nuclear matter is treated as ideal Fermi gas of the nucleons. For the 
spin- dependent operators A s the expectation value in the matter (M\A S \M) = 0, although for the 
separate polarized nucleons (N^\A S \N^) may have a nonzero value. For the operators A which do 
not depend on spin the deviation of the expectation values (M\A\M) from (0|A|0) is determined by 
incoherent sum of the contributions of the nucleons. Thus for any SU(2) symmetric spin- independent 
operator A 

(M\A\M) = (0\A\0) + p(N\A\N) (10) 
with p standing for the density of nuclear matter and 

(N\A\N) = J d 3 x(^(N\A(x)\N) - (0\A(x)\0)\ . (11) 

Since (0|A(x)|0) does not depend on x, Eq. flTTD can be presented as 

(N\A\N) = J d 3 x (N\A{x)\N) - (0\A\0) ■ V N (12) 

with Vn being the volume of the nucleon. 

The quark condensates of the same dimension d = 3 can be built by averaging of the expression 
x/jBip with B being an arbitrary 4x4 matrix over the ground state of the matter. However, any of 
such matrices can be presented as the linear combination of 5 basic matrices r a- 

Ti = I, T 2 =i fl , r 3 = 7 5 , r 4 = 7M7s, r 5 = = ^{lulu-lulu) (13) 

with I being the unit matrix. One can see, that expectation value ipT^ip vanishes in any uniform 
system, while those of ipTs^ip vanish due to conservation of parity. 

The expectation value 

*£{M\$aM M ) = Mp) (14) 

i 

takes the form v^p) = v(p)S^o in the rest frame of the matter. It can be presented as 

v(p) = E^4^' P = 5> (15) 
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with nf standing for the number of the valence quarks of the flavour "i" in the proton (neutron). 
Due to conservation of the vector current Eq.(14) presents exact dependence of this condensate on 
p. For the same reason the linear dependence on p is true in more general case of the baryon matter 

3 

Vi{p) = 2 " P> V (P) = 3 " P- ( 16 ) 

As to the expectation value (M\ipip\M), it is quite obvious that Eq. fllOD is true for the operator 
A = ipiipi if the nucleon density is small enough. The same refers to the condensates of higher 
dimension. The question is: when will the terms nonlinear in p become important ? 

Before discussing the problem we consider the lowest dimension condensates in the gas approxi- 
mation. 



2.3 Physical meaning of the scalar condensate in a hadron 



It has been suggested by Weinberg JF] that the matrix element of the operator ipiipi in a hadron is 
proportional to the total number of the quarks and antiquarks of flavour "z" in that hadron. The 
quantitative interpretation is, however, not straightforward. It was noticed by Donoghue and Nappi 
|| that such identification cannot be exact, since the operator ipiipi is not diagonal and can add 
quark-antiquark pair to the hadron. It was shown by Anselmino and Forte R ITQ] that reasonable 



assumptions on the quark distribution inside the hadron eliminate the non-diagonal matrix elements. 
However there are still problems of interpretation of the diagonal matrix elements. 

Present the quark field of any flavour 



d 3 p 



(27T) 3 (2E) 1 /2 



b s (p)u s (p)e- l ^ + df(p)v s (p)e 



-i(px) 



(17) 



with b s (p) and d+(p) eliminating quarks and creating antiquarks with spin projection s, correspond- 
ingly. This leads to 



uM) Nn P ) + v -4^N 7 (p) 



2E t {p) 



Here and N~ stand for the number of quarks and antiquarks. In the works || 10] this formula 
was analysed for the nucleon in framework of quasi-free parton model for the quark dynamics. In 
this case the normalization conditions are u s (p)u s (p) = v s {p)v s (p) = 2mj with m; standing for the 
current mass. The further analysis required further assumptions. 

In the nowadays picture of the nucleon its mass m is mostly composed of the masses of three 
valence quarks which are caused by the interactions inside the nucleon. In the orthodox nonrel- 
ativistic quark model, in which possible quark-antiquark pairs are ignored, we put Ei = rrii and 
find (N\ J^iipiipilN) = 3. In more realistic, relativistic models, there is also the contribution of the 
quark-antiquark pairs. Note also that in some approaches, say, in the bag models [11] or in the 
soliton model |I2] the motion of the valence quarks is relativistic. This reduces their contribution to 
the expectation value (N\ ^2iipiipi\N) by about 30%, since rrii/Ei < 1. 
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The conventional nowadays picture of the nucleon is that it is the system of three valence quarks 
with the constituent masses Mi ~ m/3 and the number of quark-antiquark pairs 



(N\J2AA\N) = 3 + ]T jd 3 p 



a-sip) 
2E(p) 



N s (p) 



(19) 



with a s (p) = u s (p)u s (p) = v s (p)v s (p), while N s (p) stands for the number of quark-antiquark pairs 
with momentum p. Thus, the right-hand side (rhs) of Eq.(|TU[) can be treated as the total number 
of quarks and antiquarks only under certain assumptions about the dynamics of the constituents of 
qq pairs. They should remain light and their motion should be nonrelativistic, with a s pa 2m pa 2E. 
In other models the deviation of the left-hand side (lhs) from the number 3 is a characteristic of the 
role of quark-antiquark pairs in the nucleon. 

The value of (N\ V^il-^O is related to the observables. The pion-nucleon a-term, defined by 
analogy with Eq.@ ]13| 

(20) 



1 



° = oE(^l n 5 (0)[if(0),tf(0)] \N) 



provides by using Eq.(f^) 



with 



(N\qq\N) 



2a 



m u + m d 
uu + dd . 



(21) 



qq = uu + dd . (22) 

On the other hand, |14], [15] the a-term is connected to the pion-nucleon elastic scattering amplitude 
T. Denote p, k(p', k') as momenta of the nucleon and pion before (after) scattering. Introducing 
the Mandelstam variables s = (p + k) 2 , t — (A/ — k) 2 we find the amplitude T(s,t,k 2 ,k' 2 ) in the 
unphysical point to be 

T(m 2 , 0,0; 0) = - ^ . (23) 



f 2 

The experiments provide the data on the physical amplitude 



T Una + m^) 2 , 2m 2 , m 2 , 



S 

72 

J TT 



leading to 0, |T^] 



(60 ± 7) MeV . 



(24) 



(25) 



The method of extrapolation of observable on-mass shell-amplitude to the unphysical point was 
developed by Gasser et al. [0|. They found 

a = (45 ± 7) MeV . (26) 

Note that from the point of chiral expansion, the difference E — a is of higher order, i.e. (X — a) /a ~ 

The value a = 45 MeV corresponds to (N\qq\N) pa 8. This is the strong support of the presence 
of qq pairs inside the nucleon. However direct identification of the value (N\qq\N) with the total 
number of quarks and antiquarks is possible only under the assumptions, described above. 
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2.4 Quark scalar condensate in gas approximation 



The formula for the scalar condensate in the gas approximation 

(M\qq\M) = (0\qq\0) + ^ 



m u + m d 



or 



(M\qq\M) = (0|gg|0) 1 - 



a 



P 



(27) 



(28) 



was obtained by Drukarev and Levin [^(j, |2T[. Of course, one can just substitute the semi-experimental 
value of a, given by Eq.([26|). However for the further discussion it is instructive to give a brief review 
of the calculations of the sigma-term. 

Most of the early calculations of cr-term were carried out in the framework of NJL model - 
see Eq.(l). The results were reviewed by Vogl and Weise [ 2~2| . In this approach the quarks with 
initially very small "current" masses m u « 4MeV, m d m 7MeV obtain relatively large "constituent" 
masses Mj ~ 300 — 400 MeV by four-fermion interaction, — Eq.(l). If the nucleon is treated as 
the weakly bound system of three constituent quarks, the cx-term can be calculated as the sum of 
those of three constituent quarks. The early calculations provided the value of a ~ 34 MeV, being 
somewhat smaller, than the one, determined by Eq.(^). The latter can be reproduced by assuming 
rather large content of strange quarks in the nucleon || or by inclusion of possible coupling of the 
quarks to diquarks [23j . 

In effective Lagrangian approach the Hamiltonian of the system is presented by Eq.@ with 
determined by Eq.(^) while H is written in terms of nucleon (or constituent quarks) and meson 
degrees of freedom. It was found by Gasser |24j] that 



a 



m 



dm 
dm 



with 



m, 



m u + m d 



The derivation of Eq. (|29|) is based in Feynman-Hellmann theorem |25 
Eq.(|^) is that the derivatives of the state vectors in the equation 



m = (N\H\N) 



(29) 
(30) 

The nontrivial point of 
(31) 



cancel. 



Recently Becher and Leutwyler [gfj reviewed investigations, based on pion-nucleon nonlinear 
Lagrangian. In this approach the contribution of qq pairs is 



dm dml 



m 



dml dm 



with the last factor in rhs 



dml 
dm 



ml 

7T 

m 



(32) 



(33) 
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as follows from Eq.(|8]). The calculations in this approach reproduce the value a ~ 45MeV. 



Similar calculations [27]] were carried out in framework of perturbative chiral quark model of 



Gutsche and Robson [[28[] which is based on the effective chiral Lagrangian describing quarks as 
relativistic fermions moving in effective self-consistent field. The qq pairs are contained in pions. 
The value of the a-term obtained in this model is also a ~ 45 MeV. 



30 



The Skyrme-type models provide somewhat larger values a 
. The chiral soliton model calculation gave a = 54.3 MeV pi 



50 MeV 29 and a = 59.6 MeV 



The results obtained in other approaches are more controversial. Two latest lattice QCD calcu- 
lations gave a = (18 ± 5) MeV |32| and o = (50 ± 5) MeV The attempts to extract the value 
of cr-term directly from QCD sum rules underestimate it, providing a = (25 ±15) MeV [133J and 
a = (36 ±5) MeV @. 



We shall analyse the scalar condensate beyond the gas approximation expressed by Eq. 
Subsection 2.7. 



in 



2.5 Gluon condensate 

Following Subsec.2.2 we write in the gas approximation 



CV CV C\ 

{M\—G 2 \M) = <0|— G 2 \0)+p(N\ — G 2 \N) (34) 

7T TT TT 

with notation G 2 = G^G^ . Fortunately, the expectation value (N\^-G 2 \N) can be calculated. This 
was done |5(J by averaging of the trace of QCD energy-momentum tensor, including the anomaly, 



over the nucleon state. The trace is 



= E^^-^ G2 ( 35 ) 



with b = 11 — |n, while n stands for the total number of flavours. However, (N\6£\N) does not 
depend on n due to remarkable cancellation obtained by Shifman et al. f36| 



In 

5>MMfc|JV) ~-n h (N\-^G; U G^\N) = . (36) 



h - 67V 



Here "h" denotes "heavy" quarks, i.e. the quarks, whose masses rrih are much larger than the inverse 
confinement radius /i. The accuracy of Eq.(35) is (/i/m^) 2 . Thus we only have to consider the light 
flavors u, d, s to give a reasonable approximation since m c w 1.5 GeV~ 0.3 Fm _1 . This leads to 



9 ,, T ,a s 



(N\0^\N) = -- (V|^G 2 |V) + Sm,(iV|^|iV) (37) 

P 8 7T 

with i standing for u, d and s. Since on the other hand (N\8£\N) = m one comes to 

(N\^ G 2 \N) = - ^ U-J2m(N\^i\N)\ . (38) 
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For the condensate 

g{p) = (M\ ^ G 2 \M) (39) 
Drukarev and Levin EO, pi| obtained in the gas approximation 



g(p) = g(o)-l P (m-Y,m(N\M\N)). 



(40) 



In the chiral limit m u = = and 

9{P) = g(0)-lp(m-m s (N\ss\N)y (41) 

The expectation value (iV|ss|iV) is not known definitely. Donoghue and Nappi M obtained 
(iV|ss|./V) 1 assuming, that the hyperon mass splitting in SU(3) octet, is described by the lowest 
order perturbation theory in m s . Approximately the same result (iV|ss|iV) ~ 0.8 was obtained in 
various versions of chiral perturbation theory with nonlinear Lagrangians [fTS[| . The lattice calcula- 



tions provide larger values, e.g. (N\ss\N) ~ 1.6 ||37|1 . On the contrary, the Skyrme model |12| and 



perturbative chiral quark model |27]] lead to smaller values, (iV|ss|iV) ~ 0.3. Since m ^> m s it is 
reasonable to treat the second term in the brackets of rhs of Eq. (fiTf ) as a small correction. Thus we 
can put 

g(p) = s(0) ~\pm, (42) 

which is exact in chiral SU(3) limit in gas approximation. 

One can estimate the magnitude of nonlinear contributions to the condensate g(p)- Averaging 8£ 
over the ground state of the matter one finds 

9(p) = g(0)-^(m-m s (N\ss\N))p-^e(p)p+^m s S q (p) (43) 

with e(p) standing for the binding energy of the nucleon in medium, while S q (p) denotes nonlinear 
part of the condensate (M\ss\M). One can expect the last factor to be small (otherwise we should 
accept that strange meson exchange plays large role in N — N interaction). Hence we can assume 

g(p) = g(0)-^(m + e(p))p + ^m s (N\ss\N)p (44) 

with nonlinear terms caused by the binding energy e(p). 

Thus, at least at the densities close to saturation value, corrections to the gas approximation are 
small. At p m po the value of the condensate g(p) differs from the vacuum value by about 6%. 

2.6 Analysis of more complicated condensates 

The condensates of higher dimension come from averaging of the products of larger number of 
operators of quark and (or) gluon fields. Such condensates appear also from the expansion of bilocal 
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operators of lower dimension. Say, the simplest bilocal condensate C(x) = (O| , 0(O)'0(a;)|O) is gauge- 
dependent (recall that the quarks interact with the vacuum gluon fields). To obtain the gauge- 
invariant expression one can substitute 

i/i{x) = ^(0)+x t ,D^(0) + -x fM x u D fM D l/ il;(0)+ ■■■ (45) 

with being the covariant derivatives, which replaced the usual partial derivatives [pS| . Due to 
the Lorentz invariance the expectation value C(x) depends on x 2 only. Hence, only the terms with 
even powers of x survive, providing in the chiral limit m q = 

C(x) = C(0) + x 2 ■ 1 (01-0 ^ £ GT^VIO) + ■ • • , (46) 

lO 7T 2 

where A a are Gell-Mann SU(3) basic matrices. The second term in right-hand side (rhs) of Eq . fl46|) can 
be obtained by noticing that (Ol^-D^-D^O) = ^g flu (0\ipD 2 ip\0) and by applying the QCD equation 
of motion in the form 

(&-\^ ~-mfy = 0. (47) 

The condensate (0\ip^-G^ u a^^iplO) is usually presented "in units" of (0\ipip\0), i-e. 

(0|^^GT^yV|0> = mg(0|^|0> (48) 



with mo having the dimension of the mass. The QCD sum rules analysis of Belyaev and Ioffe |39] 
gives m§ ~ 0.8 GeV 2 for u and d quarks. However instanton liquid model estimation made by Shuryak 



40[ provides about three times larger value. 

The situation with expectation values averaged over the nucleon is more complicated. There is 
infinite number of condensates of each dimension. This happens because the nonlocal condensates 
depend on two variables x 2 and (Px) with P being the four- dimensional momentum of the nucleon. 
Thus, even the lowest order term of expansion in powers of x 2 (x 2 = 0) contains infinite number of 
condensates. Say, 

(N\^(0)^(x)\N) = ^<p a ((Px),x 2 )+ix^m£ b ((Px),x 2 ) (49) 

m 

with <f(x) defined by expansion, presented by Eq.((4^). The function <p a (0, 0) is the number of the 
valence quarks of the fixed flavour in the nucleon. Presenting 

i 

$ a , b ((Px),0) = Va , b ((Px)); <Pa, b ((Px)) = I ' dae- iai * Px) <p a ^a) (50) 



we find the function <j> a {pt) to be the asymptotics of the nucleon structure function |41] and the 
expansion of (p a in powers of (Px) is expressed through expansion in the moments of the structure 
function. The next to leading order of the expansion of (p a in powers of x 2 leads to the condensate 

(^1^(0)^7,75^(0)1^) = 2P fl m^ a (51) 
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with G 



I-' 1 ' 



\£y,vafiG a a fi ■ |A a and 



£a(6) 



da9 a r b )(a,Q); 8 a{b) (a,x 2 ) 



a (b)(a,x 2 ) 
dx 2 



(52) 



The QCD sum rules analysis of Braun and Kolesnichenko [42] gave the value £ a = —0.33 GeV . 

Using QCD equations of motion we obtain relations between the moments of the functions <fi 
and 06. Denoting (F) = Jq daF(a) for any function F we find, following Drukarev and Ryskin [43 







b a a) ; (<j) b a) 



W)- i (^a) 



) = g(M 



(53) 



Situation with the nonlocal scalar condensate is somewhat simpler, since all the matrix elements 
of the odd order derivatives are proportional to the current masses of the quarks. This can be shown 
by presenting = ^(j^D + Dj^) followed by using the QCD equations of motion. Hence, in the 
chiral limit such condensates vanish for u and d quarks. The condensate containing one derivative can 
be expressed through the vector condensate and thus can be obtained beyond the gas approximation 

(M^D^M) = rmv^p) . (54) 

In the chiral limit m u = rrid = this condensate vanishes for u and d quarks. The even order 
derivatives contain the matrix elements corresponding to expansion in powers of x 2 which do not 
contain masses. In the lowest order there is the expectation value (N\ip^^G^ u a^ u ip\N) — compare 



Eq.(Hq). It was estimated by Jin et al. [44] in framework of the bag model 



(N\i>^—G<?a^\N) w 0.6 GeV 2 (55) 

7T I 



together with another condensate of the mass dimension 5 



(N\ip— — ToC^cr pvip | N) w 0.66 GeV 2 . (56) 
7T 2, 



Considering the four-quark condensates, we limit ourselves to those with colourless diquarks with 
fixed flavours. The general formula for such expectation values is 

qf = (M|^r A ^r B Vi|M) (57) 

with A,B= 1 ... 5, matrices Ta,b are introduced in Eq. (|13|) . For two lightest flavours there are thus 
5 ■ 5 • 4 = 100 condensates. Due to SU(2) symmetry Q^f = Q^f = Q AB . Due to parity conservation 
only the diagonal condensates Qfj A and also Q} 2 = Q 2 j and Q\ A = Qff have nonzero value in uniform 
matter. Since the matter is the eigenstate of the operator ^r 2 V'j, we immediately find 

Q} 2 = Pl (M\^\M) (58) 

with pi standing for the density of the quarks of i-th flavour. In the case, when the matter is composed 
of nucleons distributed with the density p, we put pi = riip with rii being the number of quarks per 
nucleon. 
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For the four-quark scalar condensate Q 11 we can try the gas approximation as the first step 
see Eq.fllO|). Using Eq.(|T2"D we find for each flavour 



(N\ipipipip\N) = I d 3 x(N\[ip(x)ip(x) - (0|^|0)] 2 |iV) + 
+ 2(Q\W\0) (N\W\N) + ^(((O^IO)) 2 - (0|#^|0>). (59) 

One can immediately estimate the second term to be about —0.09 GeV 3 . This makes the problem 
of exact vacuum expectation value to be very important. Indeed, one of the usual assumptions is 
that § 

(0|^^|0> ^ ((0|^|0)) 2 . (60) 
This means that we assume the vacuum state |0)(0| to give the leading contribution to the sum 

(0|#^|0) = 5^(0|$V|n) (n\ipip\0) (61) 

n 

over the complete set of the states |n) with the quantum numbers of vacuum. Novikov et al. K5 



showed, that Eq.(|60D becomes exact in the limit of large number of colours N c — > oo. However, the 
contribution of excited states, e.g. of the cr-meson |er)(er| can increase the rhs of Eq.(59). Assuming 
the nucleon radius to be of the order of 1 Fm we find the second and the third terms of the rhs 
of Eq.(|59D to be of the comparable magnitude. This becomes increasingly important in view of the 



only calculation of the 4-quarks condensate in the nucleon, carried out by Celenza et al. [[5J. In this 
paper the calculations in the framework of NJL model show that about 75% of the contribution of 
the second term of rhs of Eq. (BDl) is cancelled by the other ones. 



2.7 Quark scalar condensate beyond the gas approximation 

Now we denote 

(M\qq\M) = K {p) (62) 
and try to find the last term in the rhs of the equation 

2a 

«(p = k(0 + — — • p + S(p) . (63) 

m u + m d 



The first attempt was made by Drukarev and Levin [20|, ^T] in the framework of the meson-exchange 
model of nucleon-nucleon (NN) interactions. In the chiral limit m 2 — > (neglecting also the finite 
size of the nucleons) one obtains the function S(p) as the power series in Fermi momenta pp. The 
lowest order term comes from Fock one-pion exchange diagram (Fig.l). The result beyond the chiral 



limit was presented in [43 



In spite of the fact that the contribution of such mechanism to the interaction energy is a minor 
one, this contribution to the scalar condensate is quite important, since it is enhanced by the large 
factor (about 12) in the expectation value 

2m 2 

(7r|gg|7r) = 2— w 2m n ■ 12 (64) 

m u + m d 
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1 2 1 



Figure 1: The interaction of the operator qq (the dark blob) with the pion field emerging in the 
single-pion exchange. The solid lines denote the nucleons; wavy line denotes the pion. 



obtained by averaging the QCD Hamiltonian over the pion state. Using the lowest order irN coupling 
terms of the ttN Lagrangian, we obtain in the chiral limit 

S(p) = -3.2 — p (65) 
Pfo 

with pf being Fermi momentum of the nucleons, related to the density as 



2 



P = 3^ Pf ! (66) 

Pfo ~ 268 MeV is Fermi momentum at saturation point. Of course, the chiral limit makes sense only 
for p 2 F 3> m\. This puts the lower limit for the densities, when Eq.(^) is true. The value, provided 
by one-pion exchange depends on the values of ttN coupling g = gA/^fir and of the nucleon mass in 
medium. If we assume that these parameters are presented as power series in p (but not in pp) at 
low densities, the contribution of the order p 5 / 3 comes from two-pion exchange with two nucleons in 
the two-baryon intermediate state — Fig.2. 



2 

mi 



In our paper j|7 we found for p\ 3> 

S(p) = -3.2 — p - 3.1 I — j p +0(/r). (67) 




Although at saturation point m^/p F0 ~ 1/4, the discrepancy between the results of calculation of 
one-pion exchange term in the chiral limit and that with account of finite value of m\ is rather large 



43| . However, working in the chiral limit one should use rather the value of S, defined by Eq.( 



for the sigma-term, since the difference between £ and a terms contains additional powers of m, 



18] . This diminishes the difference of the two results strongly. Additional arguments in support of 



the use of the chiral limit at p close to po were given recently by Bulgac et al. pES" 



The higher order terms of the expansion, coming from the NN, NA and AA intermediate states, 
compensate the terms, presented in rhs of Eq . (|67|) to large extent. However these contributions 
are much more model-dependent. The finite size of the nucleons should be taken into account to 
regularize the logarithmic divergence. Some of the convergent terms are saturated by the pion 
momenta of the order k ~ {m{m& — m)) 1 / 2 ~ 530 MeV, corresponding to the distances of the order 
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a b c 




d e 



Figure 2: The interaction of the operator qq with the pion field created by the two-pion exchange 
between the two nucleons, denoted by the solid lines. The solid lines in the intermediate states stand 
for the nucleons or for delta isobars. The other notations are the same as in Fig.l. 




Figure 3: The behaviour of the quark scalar condensate k(p)/\k(0)\ as function of the ratio p/po 
obtained in framework of various models. The solid line shows the gas approximation |^0| , |2l| , the 
long-dashed and dashed curves present the pure p| and modified [pOl NJL model results. The 



dotted and dash-dotted lines present the result of calculation in hadronic model approach |47| . The 
dotted line corresponds to the physical value of the pion mass. The dash-dotted line shows the result 
in the chiral limit ml=0. 
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0.4 Fm, where the finite size of the nucleons should be included as well. Also, the result are sensitive 
to the density dependence of the effective nucleon mass m*. This prompts, that a more rigorous 
analysis with the proper treatment of multi-nucleon configurations and of short distance correlations 
is needed. We shall return to the problem in Sec. 4. 

Anyway, the results for the calculation of the scalar condensate with the account of the pion 
cloud, produced by one- and two-pion exchanges looks as following. At very small values of density 
P ~ Po/8, e.g. pp < only the two-pion exchanges contribute and 

S(p) = 0.8p • — . (68) 
Po 

Hence, 5* is positive for very small densities. However, for p > p we found S < 0. The numerical 
results are presented in Fig. 3. One can see the effects of interaction to slow down the tendency 
of restoration of the chiral symmetry, in any case requiring n(p) = 0. There is also the negative 
49f contribution to n(p) of the vector meson field. The sign of this term can be understood in the 



following way. It was noticed by Cohen et al. |50j that the Gasser theorem [fj4| expressed by Eq.(f2 



can be generalized for the case of the finite densities with 

S M - , (09) 

while e(p) is the binding energy. The contribution of vector mesons to rhs of Eq.(^) is 

with my standing for the vector meson mass. Since the energy caused by the vector meson exchange 

V > drops with growing my, the contribution is negative indeed. 

Another approach to calculation of the scalar condensate based on the soft pion technique was 
developed by Lyon group. Chanfray and Ericson |5l] expressed the contribution of the pion cloud 
to k(p) through the pion number excess in nuclei ||52| . The calculation of Chanfray et al. |k§ was 
based on the assumption, that GMOR relation holds in medium 

ffmf = -m (M\qq\M) . (70) 

This is true indeed, as long as the pion remains to be much lighter than the other bosonic states of 
unnatural parity. Under several assumptions on the properties of the amplitude of tcN scattering in 
medium, the authors found 

44 = - (71) 

«(0) 1 + pa/ftm* 1 ' 

and k(p) turns to zero at asymptotically large p only. This formula was obtained also by Ericson 
54| by attributing the deviations from the linear law to the distortion factor, emerging because of 



the coherent rescattering of pions by the nucleons. 

However, Birse and McGovern and Birse argued, that Eq.(|7T|) is not an exact relation 
and results from the simplified model which accounts only for nucleon-nucleon interaction, mediated 
by one pion. In framework of linear sigma model, which accounts for the irn interaction and the 
cr-meson exchange, the higher order terms of p expansion differ from those, provided by Eq. (|7T|) ■ The 
further development of calculation of the scalar condensate in the linear sigma model was made by 
Dmitrasinovic |57] . 



17 



In several works the function n(p) was obtained in framework of NJL model. In the papers of 



Bernard et al. [|5£| [59| the function /t(p) was calculated for purely quark matter. The approach was 



improved by Jaminon et al. |60| who combined the Dirac sea of quark-antiquark pairs with Fermi 
sea of nucleons. In all these papers there is a region of small values of p, where the interaction inside 
the matter is negligibly small and thus n(p) changes linearly. However, the slope is smaller than the 
one, predicted by Eq. (|2"TD . In the modified treatment Cohen et al. [|5tJ fixed the parameters of NJL 
model to reproduce the linear term. All the NJL approaches provide S > 0. 



Recently Lutz et.al |61[ suggested another hadronic model, based on chiral effective Lagrangian. 
The authors calculated the nonlinear contribution to the scalar condensate, provided by one-pion 
exchange. The value of S(po) appeared to be close to that, obtained in one-pion approximation of 



4"3"[ . Hence, all the considered hadronic models provide S < except for very small values of p. 

There is a common feature of all the described results. Near the saturation point the nonlinear 
term S(p) is much smaller than the linear contribution. Thus, Eq . (|27|) can be used for obtaining the 
numerical values of k(p) at p close to po- Hence, the condensate |/t(po)| drops by about 30% with 
respect to |«(0)|. 



3 Hadron parameters in nuclear matter 
3.1 Nuclear many-body theory 

Until mid 70-th the analysis of nuclear matter was based on nonrelativistic approach. The Schrodinger 
phenomenology for the nucleon in nuclear matter employed the Hamiltonian 

Hnr = - r^- + U(p) (72) 

and the problem was to find the realistic potential energy U(p). The deviation of nonrelativistic 
effective mass m* NR from the vacuum value m can be viewed as the dependence of potential energy 
on the value of three-dimensional momenta or "velocity dependent forces" ||62|| . The results of 
nonrelativistic approach were reviewed by Bethe |k| and by Day 



Since the pioneering paper of Walecka |JI| the nucleon in nuclear matter is treated as a relativistic 
particle, moving in superposition of vector and scalar fields V^(p) and $(p). In the rest frame of the 
matter = S^qVq and Hamiltonian of the nucleon with the three-dimensional momentum p is 

H = (ap) + (3{m + $(p)) + V (p) ■ I (73) 

with a = f° o) an d P = (o -i) being standard Dirac matrices. 

Since the scalar "a-meson" is rather an effective way to describe the system of two correlated 
pions, its mass, as well as the coupling constants of interaction between these mesons and the nucleons 
are the free parameters. They can be adjusted to fit either nuclear data or to reproduce the data on 
nucleon-nucleon scattering in vacuum. The numerous references can be found, e.g. in pBf . In both 



cases the values of Vq and $ appeared to be of the order of 300-400 MeV at the density saturation 
point. 
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The large values of the fields Vq and $ require the relativistic kinematics to be applied for the 
description of the motion of the nucleons. 

In the nonrelativistic limit the Hamiltonian (|73|) takes the form of Eq.(|72|) with m* NR being 
replaced by Dirac effective mass m*, defined as 

m* = m + $ (74) 

and 

U = V + $ . (75) 

At the saturation point the fields Vq and $ compensate each other to large extent, providing U ~ 
—60 MeV. This explains the relative success of Schrodinger phenomenology. However, as shown by 
Brockmann and Weise H67H , the quantitative description of the large magnitude of spin-orbit forces 
in finite nuclei requires rather large values of both $ and Vq. 

In the meson exchange picture the scalar and vector fields originate from the meson exchange 
between the nucleons of the matter. The model is known as quantum hadrodynamics — QHD. In 
the simplest version (QHD-1) only scalar cx-mesons and vector u mesons are involved. In somewhat 
more complicated version, known as QHD-2 [B8| some other mesons, e.g. the pions, are included. 
The matching of QHD-2 Lagrangian with low energy effective Lagrangian was done by Furnstahl 
and Serot |69|| . 

The vector and scalar fields, generated by nucleons, depend on density in different ways. For the 
vector field 

V(p) = N v (p)g v 6( PF -p) (76) 

with gv the coupling constant, while Ny = (unIqu^) /2E with un(p) standing for nucleon bispinors 
while 

s(p,p) = V (p) + (p 2 + m* 2 (p)) 1/2 . (77) 

One finds immediately that Ny = 1 and thus V(p) is exactly proportional to the density p. On the 
other hand, in the expression for the scalar field 

Hp) = 4 J -0^ N s (p) ■ g s 6( PF - p) (78) 

the factor iV s = m* / E. Thus, the scalar field is a complicated function of density p. 

The saturation value of density po can be found by minimization of the energy functional 

1 " 



£{p) = - / e F (p)dp (79) 







with e F (p) = z(pf,p) being the single-particle energy at the Fermi surface. Thus, in QHD the 
saturation is caused by nonlinear dependence of the scalar field $ on density. 

The understanding of behaviour of axial coupling constant in nuclear matter (?a(p) requires ex- 
plicit introduction of pionic degrees of freedom. The quenching of gA at finite densities was predicted 
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by Ericson ||7(J from the analysis of the dispersion relations for irN scattering. The result was con- 
firmed by the analysis of experimental data on Gamow-Teller /3-decay of a number of nuclei carried 
out by Wilkinson |7l] and by investigation of beta decay of heavier nuclei — see, e.g., fF2 



g A (0) = 1.25 ; g A {p ) = 1.0 . (80) 

The quenching of g A as the result of polarization of medium by the pions was considered by Ericson 
et al. f73|. The crucial role of isobar-hole excitations in this phenomena was described by Rho [74 



Turning to the characteristics of the pions, one can introduce effective pion mass m* by considering 
the dispersion equation for the pion in nuclear matter (see, e.g., the book of Ericson and Weise |75| ): 

to 2 - k 2 - U p (lo, k)-m* 2 = . (81) 

Here ui and k are the pion energy and three-dimensional momenta, n p is the p-wave part of the pion 
polarization operator. Hence, U p contains the factor k 2 . The pion effective mass is 

mf = ml + U s (u,k) (82) 

with n s being the s-wave part of polarization operator. 

Polarization operator n s (as well as U p ) is influenced strongly by the nucleon interactions at 
the distances, which are much smaller than the average inter-nucleon distances ~ Tn~ . Strictly 
speaking, here one should consider the nucleon as a composite particle. However, there is a possibility 
to consider such correlations in framework of hadron picture of strong interactions by using Finite 
Fermi System Theory (FFST), introduced by Migdal |[76|| . In framework of FFST the amplitudes of 
short-range baryon (nucleons and isobars) interactions are replaced by certain constant parameters. 
Hence, behaviour of m* can be described in terms of QHD and FFST approaches. 

As well as any model based on conception of NN interaction, QHD faces difficulties at small 
distances. The weak points of the approach were reviewed by Negele [[77j and by Sliv et al. J78| 



Account of the composite structure of nucleon leads to the change of some qualitative results. Say, 
basing on the straightforward treatment of the Dirac Hamiltonian Brown et al. |79| found a significant 



term in the equation of state, arising from virtual NN pairs, generated by vector fields. The term 



would have been important for saturation. However, Jaroszewicz and Brodsky |5(J and also Cohen 
8l| found that the composite nature of nucleon suppresses such contributions. 



Anyway, to obtain the complete description, we need a complementary approach, accounting for 
the composite structure of hadrons. For pions it is reasonable to try NJL model. 



3.2 Calculations in Nambu— Jona— Lasinio model 

In NJL model the pion is the Goldstone meson, corresponding to the breaking of the chiral symmetry. 
The pion can be viewed as the solution of Bethe-Salpeter equation in the pseudoscalar quark- 
antiquark channel. The pion properties at finite density were investigated in frameworks of SU(2) 
and SU(3) flavour NJL model ]58|, |59[| . It was found that the pion mass m* (p) is practically constant 
at p < po, increasing rapidly at larger densities, while f*(p) drops rapidly. These results were obtained 
rather for the quark matter. Anyway, as we mentioned in Subsec.2.7, at small p the condensate K.(p), 
obtained in this approach, does not satisfy the limiting law, presented by Eq. (|27|) . 
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However, the qualitatively similar results were obtained in another NJL analysis, carried out by 
Lutz et al. fl82| . The slope of the function n{p) satisfied Eg. fl27|) . The pion mass m*(p) increased 



with p slowly, while /*(p) dropped rapidly. The in- medium GMOR relation, expressed by Eq.( 
was satisfied as well. 



Jaminon and Ripka [83] considered the modified version of NJL model, which includes the dilaton 
fields. This is the way to include effectively the gluon degrees of freedom. The results appeared to 
depend qualitatively on the way, the dilation fields are included into the Lagrangian. The pion mass 
can either increase or drop with growing density. Also the value of the slope of n(p) differs strongly 
in different versions of the approach. In the version, which is consistent with Eq.(|2TD the behaviour 
of /*(p) and m^(p) is similar to the one, obtained in the other papers, mentioned in this subsection. 

Note, however, that the results which predict the fast drop of f*(p) have, at best, a limited region 
of validity. This is because the pion charge radius r n is connected to the pion decay constant by the 



relation obtained by Carlitz and Creamer [|4 



<^> 1/2 = ^ ^ 

providing (r 2 ) 1 / 2 w 0.6 Fm. Identifying the size of the pion with its charge radius, we find that at 
(r 2 (p)) 1//2 becoming of the order of the confinement radius r c ~ 1 Fm, the confinement forces should 
be included and straightforward using of NJL is not possible any more. Thus, NJL is definitely not 
true for the densities, when the ratio f*(p)/ f n becomes too small. Anyway, one needs 

^ > 0.6 . (84) 



For the results, obtained in [82] this means that they can be true for p < 1.3po only. 



3.3 Quark— meson models 



This class of models, reviewed by Thomas |TT| is the result of development of MIT bag model, 
considering the nucleon as the system of three quarks in a potential well. One of the weak points of 
the bag-model approach is the absence of long-ranged forces in NN interactions. In the chiral bag 
model (CBM) the long-ranged tail is caused by the pions which are introduced into the model by 
requirement of chiral invariance. In the framework of CBM the pions are as fundamental degrees of 
freedom as quarks. In the cloudy bag model these pions are considered as the bound states of qq 
pairs. The model succeeded in describing the static properties of free nucleons. 

Another model, suggested by Guichon |£J is a more straightforward hybrid of QHD and QCD. 
The nucleon is considered as a three-quark system in a bag. The quarks are coupled to a- and io- 
mesons directly. Although this quark-meson coupling model (QMC) was proposed by its author as "a 
caricature of nuclear matter" , it was widely used afterwards. The parameters of a- and cu-mesons and 
the bag radius, which are the free parameters of the model were adjusted to describe the saturation 
parameters of the matter. The fields $ and V appear to be somewhat smaller than in QHD. Thus, 
the values of m*/m and g\/g are quenched less than in QHD [j36 |. On the other hand, the unwanted 



iVW pairs are suppressed. The nonlinearity of the scalar field is the source of saturation. 
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The common weak point of these models are well known. Say, there is no consistent procedure to 
describe the overlapping of the bags. It is also unclear, how to make their Lorentz transformations. 



3.4 Skyrmion models 

This is the class of models with much better theoretical foundation. They originate from the old 
model, suggested by Skyrme |$7fl. The model included the pions only, and the nucleon was the 



soliton. Later Wess and Zumino [38] added the specific term to the Lagrangian, which provided the 
current with the non-vanishing integral of the three-dimensional divergence. That was the way, how 
the baryon charge manifested itself. 

Thus, in framework of the approach most of the nucleon characteristics are determined by Dirac 
sea of quarks and by the quark-antiquark pairs, which are coupled into the pions. The model can 
be viewed as the limiting case R — > of the chiral bag model, where the description in terms of the 
mesons at r > R is replaced by description in terms of the quarks at r < R f89| . 

In the framework of the Skyrme model Adkins et al. |]90| , calculated the static characteristics 
of isolated nucleons. A little later Jackson et al. [|9l| investigated NN interaction in this model. 



The model did not reproduce the attraction in NN potential. It was included into modified Skyrme 



Lagrangian by Rakhimov et al. [93 in order to calculate the renormalization of gA, m and f n in 



nuclear matter. The magnitude of renormalization appeared to be somewhat smaller than in QHD. 



The approach was improved by Diakonov and Petrov — see a review paper |j4j] and references 
therein. The authors build the chiral quark-soliton model of the nucleon. It is based on quark- 
pion Lagrangian with the Wess-Zumino term and with spontaneous chiral symmetry breaking. The 
nucleon appeared to be a system of three quarks, moving in a classical self-consistent pion field. The 
approach succeeded in describing the static characteristics of nucleon. It provided the proper results 
for the parton distributions as well. However, the application of the approach to description of the 
values of nucleon parameters in medium is still ahead. 

3.5 Brown-Rho scaling 



Brown and Rho [951 assumed that all the hadron characteristics, which have the dimension of the 



mass change in medium in the same manner. The universal scale was assumed to be 

X (p) = (-n(p)) 1/3 . (85) 
Thus, the scaling which we refer to as BR1 is 



m U x(0) 

The pion mass was assumed to be an exception, scaling as 

™*M = fx(p)\ 1/2 
m * \x(o)y 



6) 



(87) 
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Thus, BR1 is consistent with in-medium GMOR relation. Also, in contrast to NJL, the pion mass 
drops with density. 

Another point of BR1 scaling is the behaviour 

9*a(p) = 9a{0) = const . (88) 



Consistency of Eqs. (|80D and (88) can be explained in such a way. Renormalization expressed by 



Eq.(jSOf) is due to A-hole polarization of medium. It takes place at moderate distances of the order 
m~ , reflecting rather the properties of the medium, but not the intrinsic properties of the nucleon, 
which are discussed here. 



Another version of Brown-Rho scaling [Q, which we call BR2 is based on the in-medium GMOR 
relation, expressed by Eq. (|70|) . It is still assumed that 

^ = f, (89) 

but the pion mass is assumed to be constant 

m* « m, , (90) 



(91) 



and thus 

K (<p±) 1/2 
U U(o)/ 

instead of 1/3 law in BR1 version — Eq . (|3~6l) . Note, however, that assuming m*(po) = l-OSm^ |)6| 
we find, using the results of subsection 2.6 

= 0.76 . (92) 

Jtt 



This is not far the limit determined by Eq. (|8"3|) . At larger densities the size of the pion becomes of 
the order of the confinement radius. Here the pion does not exist as a Goldstone boson any more. In 
any case, some new physics should be included at larger densities. If Eq.(^) is assumed to be true, 
this happens at p « 1.6po- 



3.6 QCD sum rules 

In this approach we hope to establish some general relations between the in- medium values of QCD 
condensates and the characteristics of nucleons. 

The QCD sum rules were invented by Shifman et al. || and applied for the description of the 
mesonic properties in vacuum. Later Ioffe |57J expanded the method for the description of the 
characteristics of nucleons in vacuum. The main idea is to build the function G(q 2 ) which describes 
the propagation of the system ("current") with the quantum numbers of the proton. (The usual 
notation is II(g 2 ). We used another one to avoid confusion with pion polarization operator, expressed 
by Eq.flSTD). The dispersion relation 
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is considered at q 2 — > — oo. Imaginary part in the rhs is expressed through parameters of observable 
hadrons. Due to asymptotic freedom of QCD lhs of Eq. (|93|) can be presented as perturbative series 
in —q 2 with QCD vacuum condensates as coefficients of the expansion. Convergence of the series 
means that the condensates of lower dimension are the most important ones. 

The method was used for the calculation of characteristics of the lowest lying hadron states. 
This is why the "pole+continuum" model was employed for the description of lmG(k 2 ) in the rhs of 
Eq. (^) . This means that the contribution of the lowest lying hadron was treated explicitly, while all 
the other excitations were approximated by continuum. In order to emphasise the contribution of the 
pole inverse Laplace (Borel) transform was applied to both sides of Eq. fl93| ) in the papers mentioned 
above. The Borel transform also removes the polynomial divergent terms. 



Using QCD sum rules Ioffe |97| found that the nucleon mass vanishes if the scalar condensate 
turns to zero. Numerically, p9| , |9~T| , |98 



rn 



-2(2vr) 2 (0|gg|0) 



1/3 



(94) 



Later the method was applied by Drukarev and Levin |20| , [21], [99|] for investigation of properties of 
nucleons in the nuclear matter. The idea was to express the change of nucleon characteristics through 
the in-medium change of the values of QCD condensates. The generalization for the case of finite 
densities was not straightforward. Since the Lorentz invariance is lost, the function G m (q) describing 
the propagation of the system in medium depends on two variables, e.g. G m = G(q 2 , q Q ). Thus, each 
term of expansion of G m in powers of q~~ 2 may contain infinite number of local condensates. In the 
rhs of dispersion relation it is necessary to separate the singularities, connected with the nucleon 
from those, connected with excitation of the matter itself. 

We shall return to these points in Sec. 5. Here we present the main results. The method provided 
the result for the shift of the position of the nucleon pole. The new value is expressed as a linear 
combination of several condensates with vector condensate v(p) and scalar condensate n(p) being 



most important 20, 21, 9£ 



m m = m + CiK,(p) + C 2 v(p) . (95) 

On the other hand, 

m m -m = u(l + o(—)) (96) 



with U being single-particle potential energy of the nucleon. Hence, the scalar forces are to large 
extent determined by the cx-term. 

The Dirac effective mass was found to be proportional to the scalar condensate 

m*(p) = K(p)F(p) (97) 

with F(p) containing the dependence on the other condensates, e.g. on vector condensate v(p). 

Using Eqs.(|T5D and ( |i~6"| ) we see, that v(p) is linear in p. Thus, the main nonlinear contributions to 
the energy S(p) presented by Eq.([79f) come from nonlinearities in the function n(p). For the saturation 
properties of the matter the sign of the contribution S(p) becomes important. The nonlinearities of 
the condensate n{p) can be responsible for the saturation if S < 0. Calculations of Drukarev and 



Ryskin |43[ show that the saturation can be obtained at reasonable values of density with reasonable 
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value of the binding energy. Of course, this result should not be taken too seriously, since it is very 
sensitive to the exact values of cr-term. It can be altered also by the account of higher order terms. 
(However, as noted by Birse [56|, the QHD saturation picture is also very sensitive to the values of 
the parameters 
Lutz et.al 161. 



. Similar saturation mechanism was obtained recently in the approach developed by 
Anyway, it can be a good starting point to analyse the problem. 



4 First step to self-consistent treatment 

As we have seen in Sec. 2 in the gas approximation the scalar condensate n(p) is expressed through 
the observables. However, beyond the gas approximation it depends on a set of other parameters. 
Here we show how such dependence manifests itself in a more rigorous treatment of the hadronic 
presentation of nuclear matter. 



4.1 Account of multi-nucleon effects in the quark scalar condensate 

Now we present the main equations, which describe the contribution of the pion cloud to the con- 
densate k(p). Recall, that the pions are expected to give the leading contribution to the nonlinear 
part S(p) due to the large expectation value (vr|gg|7r) — Eq.( 



In order to calculate the contribution we employ the quasiparticle theory, developed by Migdal 
for the propagation of pions in matter [ |100| |. Using Eq. (|5UD , we present S(p) through the derivative 



of the nucleon self-energy with respect to m 



2 . 



B 



Sb = -C B T I ( g 3 (r|D V, k)g B (p — k) — r»Ug( Wj k)g%(p - k)) . (98) 

Here B labels the excited baryon states with propagators g B and ttNB vertices T B . The pion 
propagator D includes the multi-nucleon effects (D~ l is the lhs of Eq.([H])). The second term of the 
rhs of Eq. (E8) , with the index " 0" corresponding to the vacuum values, subtracts the terms, which 
are included into the expectation value already. The coefficient C B comes from summation over the 
spin and the isospin variables. Integration over nucleon momenta p is limited by the condition p < pp. 
The factor T stands for the expectation value of the operator qq in pion, i.e. T = (7r|gg|7r) = m^/rh. 
Of course, Eq . fl98|) , illustrated by Fig. 4 corresponds to the Lagrangian which includes the lowest 
order irN interactions only. 



The pion propagator in medium can be viewed as the solution of the Dyson equation ]75, [76[] 
Fig.5: 

D = D + D TID (99) 

with 

1%, k) = 4tt J 0^ A(p>; cu, k) , (100) 
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1 2 1 

a 




d e 



Figure 4: a)The interaction of the operator qq (the dark blob) with the pion field. The solid 
line denotes the nucleon; the wavy line stands for the pion; b,c) The diagrammatic presentation of 
Eq.(98) with the nucleon in the intermediate state. The bold wavy line denotes the pion propaga- 
tor renormalized due to baryon-hole excitations in the framework of FFST; d,e)The diagrammatic 
presentation of Eq.(98) with the A-isobar (double solid line)in the intermediate state. 



VWW =wwv + vwQwv 



Figure 5: The Dyson equation (99) for the pion propagator in medium in the quasiparticle-hole 
formalism. Wavy line denotes the vacuum pion propagator, bold wavy line stands for the propagator 
in matter. The dark angle denotes the correlations. 
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while A(p'; u, k) stands for the amplitude of the forward ttN scattering (all the summation over the 
spin and isospin variables is assumed to be carried out) on the nucleon of the matter with the three- 
dimensional momentum p' . Of course, the pion is not on the mass-shell. 

Neglecting the interactions inside the bubbles of Fig. 5 (this is denoted by the upper index "(0)") 
we can present 

a<°> = E4 0) ; n(0) = E4 0) 

A° B = c B f 2 B (k)A B {p';u;,k) (101) 

with cb being a numerical coefficient, 

A B (p';uj,k) = g B (e' + oj,p' + k)+g B (e' -oj,p' -k) . (102) 

The factors T B (k) come from the vertex functions. Considering p-wave part of polarization operator 
only (the s-wave part is expressed through the pion effective mass m* — Eq. fl82D ), we present 

r|(*0 = ~9l NB k 2 d 2 NB (k) (103) 

with d 2 NB accounting for the finite size of the baryons, g n NB is the coupling constant. 

Starting the analysis with the contribution of the nucleon intermediate state (B = N) to Eqs. fllOll ) 
and ( |102|) , we see that in the nonrelativistic limit we can present the first term in rhs of Eq.( |102| ) as 



iKV + aJ + k) = ^ + ^-^> (104) 

LO + E p i — 6 p i + k 

(similar presentation can be written for the second term) with e q = q 2 /2m* + U, while U stands for 
the potential energy. Hence, the terms, containing U cancel and all the dependence on the properties 
of the matter enters through the effective mass m*. This enables to obtain the contribution to the 
polarization operator 

n!? = -49l NN k 2 d% N (k) ^(^k) ( 105 ) 

with explicit analytical expression for <p$(u,k) presented in [75, 104 1 , the static long-wave limit is 
0^(0,0) = 1. 

Such approach does not include the particle-hole interactions in the bubble diagram of Fig. 5. 
The short-range correlations can be described with the help of effective FFST constants, as it was 
mentioned above. Using the Dyson equation for the short-range amplitude of nucleon-hole scattering 
one finds 

ITv = -4gl NN k 2 d% N (k) M",k) (106) 

with 

AO) 



= : r \L . ( 107 ) 



if only the nucleon intermediate states are included. 

The long-ranged correlations inside the bubbles were analysed by Dickhoff et al. [|101|| . It was 
shown, that exchange by the renormalized pions inside the bubbles ("bubbles in bubbles") can be 
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accounted for by the altering of the values of FFST constants. The change in the numerical values 
does not appear to be large. 



The usual approach includes also the A-isobar states in the sums in Eq. (|98|) . Until the particle- 
hole correlations are included, the total p-wave operator IT( ) is just the additive sum of the nucleon 
and isobar terms, i.e. II^ ) = ITjy + n^. Also, one can obtain analytical expression, similar to 
Eq. ( |105| ) for the contribution under a reasonable assumption on the propagation of A-isobar in 
medium (see below). However, account of the short-range correlations makes the expression for the 
total p-w&ve polarization operator more complicated. We use the explicit form presented by Dickhoff 
et al. |0| 

n = nTv + riA 



with 



11* = 11$ [ 1 - ( 7 a - 7aa)^- 1 /E (108) 

n (o) ' 

n A = ni 0) ( 1 + ( 7 A - Inn)-^- ) /E . (109) 



Denominator E has the form 



n (o) n (o) n (o) n (o) 

E = 1 - lNN L -f- - 7aa^- + (7^7aa - 7a) • (HO) 

The effective constants T are related to FFST parameters g' as follows: 

n 9nn . „, n 9na „, n 9aa fiii\ 

7NN = ^0^2 ' TA = Cq- Z , 7AA = > (HI) 

9nNN g-nNNg-KNA 9-kNA 

where Co is the normalization factor for the effective particle-hole interaction in nuclear matter. We 
use Co = n 2 /pFm*, following J75J. (Note, that there is some discrepancy in the notations used by 
different authors. Our parameters 7 coincide with those, used in [|102j| . We use the original FFST 
parameters g' of J75J, which are related to the constants G' Q of [102] as g' = G'/2). The short-range 



interactions require also renormalization of the vertices T 2 NB — > T 2 NB x 2 NB with 

n (o)\ 1 / jjCo) 1 

XnNN = ( 1 + (7a - 7aa)-^- I / E ; x nNA = I 1 + (ta - Jnn)-^- ] / E . (112) 



In our paper [|103|1 we calculated the contribution S(p), presented by Eq.(^), using nucleons and 
A-isobars as intermediate states. The integration over uj requires investigation of the solutions of the 
pion dispersion equation-Eq. (|8l|) . 



4.2 Interpretation of the pion condensate 



The pion dispersion equation [75, 101 



is 



m* n 2 + k 2 [l + x (u,k) 



(113) 
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with the function x introduced as U p (u),k) = —k 2 x{u,k). It is known to have three branches of 
solutions U)i(k) (classified by the behaviour of the functions Wi(k) at k — > 0). If the function x( U) 7 ^) 
includes nucleons only as intermediate states and does not include correlations, we find k 2 \ — > 
at k — > 0. This is the pion branch for which clv(0) = m*. If the correlations are included, the 
denominator in the rhs of Eq. (|107|) may turn to zero at k — > 0, providing the sound branch with 
u) s (0) = 0. Inclusion of A-isobars causes the contribution to x{ u > proportional to [m^ — m — uj]" 1 . 
Thus, there is a solution with o>a(0) = — m, called the isobar branch. 

The trajectories of the solutions of Eq.(|8lD on the physical sheet of Riemann surface were studied 



by Migdal [ 104 |. Their behaviour on the unphysical sheets was investigated recently by Sadovnikova 
[105 1 and by Sadovnikova and Ryskin [|106fl . In these papers it was shown, that besides the branches, 
mentioned above, there is one more branch starting from the value o> c (0) = m* and moving on the 
unphysical sheet for larger k > 0. The branch comes to the physical sheet at certain value of k if the 
density exceeds certain critical value pc- Here u c is either zero or purely imaginary and thus u 2 < 0. 
(However, this is true if the isobar width = 0, for the finite values of Ta we find u c to be complex 



and Reu 2 < 0.) This corresponds to the instability of the system first found by Migdal [ 100|| and 



called the "pion condensation". On the physical sheet u c {k) coincides with the solutions, obtained 
10011 , [|104| . However, contrary to ||100|| , ||104|| , the u c (k) is not the part of zero-sound branch. 

To follow the solution uj c (k), let us present the function which enters Eq.( |105| ) as 



in 



^\u,k) = ^\u,k) + ^\-u,k) (114) 
with the explicit expression for < k < 2pp 



-urn* + kp F (kp F ) 2 - (uom* - k 2 /2)' 




^ * ln( (n5) 

1 1 ' ' " 1 1 \um-kpF + k 2 /2 J J K J 



At k > 2pp the expression for ip N (oo, k) takes another form (see ||104|| ) but we shall not need it here. 



It was shown in ||105|| , ||106|| that, if the density p is large enough (p > p c ), there is a branch of 



solutions uj 2 (k) < 0, which is on the physical sheet for certain interval k\ < k < k 2 of the values of 
k. At smaller values k < k\ the branch goes to the unphysical sheet through the cut 



< u < — [PF-- ), (H6) 




generated by the third term in the rhs of Eq. (|115|) . At larger values of k > k 2 the solution u c goes 



away to the unphysical sheet through the same cut. The zero-sound wave goes to the unphysical 
sheet through another cut: 

k ( k\ k ( k\ . . 

-,\vr-i)< - < d") 

caused by the second term in the rhs of Eq. ( |115| ) . 

The value of the density pc, for which the solution ui c penetrates to the physical sheet, depends 
strongly on the model assumptions. Say, if the contribution of isobar intermediate states is ignored, 
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the value of pc is shifted to unrealistically large values pc > 25po- Inclusion of both nucleon and 
isobar states and employing of realistic values of FFST constants leads to pc ~ l-4po under additional 
assumption m*(p) = m n (0). 

The zero values of u c (k) at certain nonzero values of k signals on the instability of the ground 
state . New components, like baryon-hole excitations with the pion quantum numbers emerge in the 
ground state of nuclear matter. Thus, the appearance of the singularity corresponding to u 2 = 
shows, that the phase transition takes place. 

Note, however, that the imaginary part of the solution u c {k) is negative. Thus, there is no 
" accumulation of pions" in the symmetric nuclear matter, contrary to the naive interpretation of the 
pion condensation. 

The situation is much more complicated in the case of asymmetric nuclear matter. In the neutron 
matter the instability of the system emerges at finite values of u, because of the conversion n — > 
p + 7r~ jnj. This process leads to the real accumulation of pions in the ground state. In the charged 
matter with the non-zero value of the difference between the neutron and proton densities there is 
an interplay of the reactions n p + tt~ and beta decays of nucleons. 



4.3 Quark scalar condensate in the presence of the pion condensate 



Now we turn back to the calculation of the condensate k(p). Note first, that if the isobar width is 
neglected, we find n(p) — > +00 at p — ► pc- The reason is trivial. When p — > pc, the contribution 
S(p), described by Eq.(OH) becomes 



S 



dud 3 k 



[uj 2 -u 2 (p,k)\ 2 

The curve uj c (p c , k) turns to zero at certain k = k c , being uo c = a(k — k c ) 2 at \k — k c \ <C k c . Thus 

duk 2 d\k — k c \ 



(118) 



S 



[u 2 - a 2 (k - k c ) 2 f 



00 . 



(119) 



Hence, S(pc) = +00 and n(pc) = +00. Once /c(0) < 0, we find that at certain p c h < pc the 
scalar condensate n(p) turns to zero. This means that the chiral phase transition takes place before 
the pion condensation. (We shall not discuss more complicated models, for which the condition 
K{pch) = is not sufficient for the chiral symmetry restoration). At larger densities the pion does 
not exist any more as a collective Goldstone degree of freedom. Also the baryon mass vanishes (if 
very small current quark masses is neglected), and we have to stop our calculations, based on the 
selected set of Feynman diagrams (Fig.4) with the exact pion propagator. 

The in-medium width of delta isobar Ta (the probability of the decay) depends strongly on the 
kinematics of the process. We can put T&(uj,k) = due to the limitation on the phase space of 
the possible decay process 



103 



106] 



we 



Thus, following the paper of Sadovnikova and Ryskin 
find that the chiral symmetry restoration takes place at the densities, which are smaller, than those 
corresponding to the pion condensation: 



Pch < Pc ■ 



(120) 
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This means, that the pion condensation point cannot be reached in framework of the models, which 
do not describe the physics after the restoration of the chiral symmetry. 



4.4 Calculation of the scalar condensate 
4.4.1 Parameters of the model 

Now we must specify the functional dependence and the values of the parameters which are involved 
into the calculations. The ttNN coupling constant is 

9 ^ NN ~2m~ = 2U ' (121) 

- see Eq.(f|). The 7riVA coupling constant is 

9nNA = Ca9ttNN (122) 

with the experiments providing c A ~ 2 |75fl . This is supported by the value c A ~ 1.7, calculated in 
the framework of Additive Quark Model (AQM). 

The form factor djvs(^) which enters Eq. (|103|) is taken in a simple pole form [75 



1 - ml/ A 



dNB ~ 1 + k 2 /A 2 B (123) 
with A N = 0.67 GeV, A A = 1.0 GeV. 

We use mostly the values of FFSI parameters, presented in |J107|j: g' NN = 1.0, g' NA = 0.2, g' AA = 0.8 



- referring to these values as to set "a". We shall also check the sensitivity of the results to the 
variation of these parameters. 

It is know from QHD approach that the nucleon effective mass may drop with density very rapidly. 
Thus, we must adjust our equations for description of the case, when the relativistic kinematics 
should be employed. We still include only the positive energy of the nucleon propagator, presented 
by Eq.( |104|) . However we use the relativistic expression for 



e p -e p+k = ^p 2 + m* 2 - ^(p + k) 2 + m* 2 . (124) 

The propagator of A-isobar is modified in the same way. The explicit equations for the functions 
$at Aj accounting for the relativistic kinematics are presented in ||103| . 



4.4.2 Fixing the dependence m*(p) 

As we have seen above, the contribution of nucleon-hole excitations to S(p) depends explicitly on the 
nucleon effective mass m*(p). Here we shall try the models, used in nuclear physics, which determine 
the direct dependence m*(p). One of them is the Fermi liquid model with the effective mass described 
by Landau formula [M , PI , jT09 



m*(p) _ i/ (i + 2mpF 



m /v 7r 2 



/i) ■ (125) 
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In QHD approach the effective mass m* is the solution of the equation pS 



m = m — cm 



( 2 , *2^1/2 ,2, PF + (P% + m* 2 ) 1/2 

Pf{Pf + m ) ~ m hi 



mr 



(126) 



corresponding to the behaviour 

m* = m (l - f 2 p) (127) 

in the lowest order of expansion in powers of Fermi momentum pp. The coefficients in Eqs. (|125|) , 
(|127|) can be determined by fixing the value m*(p ). 

Assuming, that all the other parameters are not altered in medium: /* = f n , m A — m* = m A — m, 
m* = uin, c A = ca, g\ = 9a(po) ~ 1-0, we find the point of chiral symmetry restoration to depend 
strongly on the value m*(po), being stable enough under the variation of FFST parameters and of 
the parameter ca — Fig. 6. Fixing m*(po) = 0.8m, we find p c h < po, in contradiction to experimental 
data. Even in a simplified model with the width of A-isobar being accepted to coincide with its 
vacuum value Ta = 115 MeV, we find p c h ~ 1.15po- The value |re(po)| <S l K (0)| looks unrealistic, since 
there are practically no strong unambiguous signals on partial restoration of the chiral symmetry 



at the saturation value of density po |56| . Hence, here we also come to contradiction with the 
experimental data. 

The situation is less critical for the smaller values of m*(po) /m. Say, for Ta = 115 MeV we find 
p c h = 1.7 po, assuming m*(p )/m = 0.7. However, under realistic assumption Ta = we come to 
Pch < Po- 

Note, that there is another reason for the point of the pion condensation to be unaccessible by 
our approach. The perturbative treatment of ttN interaction becomes invalid for large pion fields. 
In the chiral ttN Lagrangians the nN interaction is described by the terms of the type 

L nN = $U+(i lfl d»)U*l; (128) 

with 

i 

U = exp— 7 5 (r<p) . 

The conventional version of pseudovector ttNN Lagrangian employed above may be treated either 
as the lowest term of expansion of the matrix U in powers of the ratio <p//?r (identifying the pion with 
<p-field) or as the interaction with the field (p = /„- sin ((r<p) //„■). In any case, the whole approach is 
valid only, when the pion field is not too strong (<p < f m or (p < f n correspondingly). However, the 
strict quantitative criteria for the region of validity of Eq.(p8|) is still obscure. 

The strong dependence of the results on the value of m*(po)/m forces us to turn to self-consistent 
treatment of the hadron parameters and the quark condensates. 



4.4.3 Self-consistent treatment of nucleon mass and the condensate 

Now we shall carry out the calculations in framework of the model, where the nucleon parameters 
depend on the values of condensates. In other words, instead of the attempt to calculate the con- 
densate k(p, Ui{p)) with Hi standing for the hadron parameters (yi = m* N , m* A , /*, . . .), we shall try 



32 




Figure 6: The function k(p)/\k(0)\. The solid curve presents the result obtained with g' NN = 1.0, 
9na = 0-2) A'aa = 0-8, c a = 2.0, Ta = 0.115 GeV and nucleon effective mass given by Eq.( |125| ) 
The other curves are obtained with the values of some of the parameters or the shape of the density 
dependence of the effective mass being modified: a) Dependence of k(p)/\k(0)\ on the variation of 
nuclear parameters. The dashed curve corresponds to the calculation with ca = 1.7, the dotted curve 
~ t° 5 , aa = 1-2 ) dot-dashed curve - to g' NN = 0.7. b) Dependence of n(p)/\n(0)\ on the isobar width. 
Dotted curve corresponds to the calculation with Ta = 0.07 GeV, dot-dashed curve - to Ta = 0.05 
GeV, dashed curve - to Ta = 0.01 GeV. c) Dependence of n(p)/\n(0)\ on the shape of m*(p). 
Dashed curve corresponds to Walecka formula ( |127| ) with m* (m*(p = pa) = 0.8m). Dot-dashed 
curve is obtained in framework of Walecka model with m*(p = p ) = 0.7m. 
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to solve the equation 

k( P ) = K(p,yMp)MP))) ( 129 ) 

with Cj(p) standing for the other QCD condensates. Here K is the rhs of Eq.(62). Strictly speaking, 
we should try to obtain similar equations for the condensates Cj(p). 

We shall assume the physics of nuclear matter to be determined by the condensates of lowest 
dimension. In other words, we expect that only the condensates, containing the minimal powers of 
quark and gluon fields are important. The condensates of the lowest dimension are the vector and 



scalar condensates, determined by Eqs. (|Tq) and (p2j) and also the gluon condensate — Eq . (p4|) . As 



we saw in Subsect.2.7, the relative change of the gluon condensate in matter is much smaller than 
that of the quark scalar condensate. Thus we assume it to play a minor role. Hence, the in- medium 
values of n(p) and v(p) will be most important for us, and we must solve the set of equations 

K(p) = K{p,yi) 

V% = Vi(v(p),K(p)) ■ (130) 

Fortunately, the vector condensate v(p) is expressed by simple formulas fll^D and (|IED due to the 
baryon current conservation. 

The idea of self-consistent treatment is not a new one. Indeed, Eqs. (2) and (3) provide an 
example of Eq. ( |130| ) for N JL model in vacuum, with the only parameter t/i = m. 

As to parameters y^ which are m* ,m* A , f*, etc., there are several relations which are, to large 
extent, model-independent. Besides the in-medium GMOR relation — Eq . (|70|) , we can present in- 
medium GT relation 

9*nNN = 7fJ~ ■ ( 131 ) 

Recalling that GT relation means, that the neutron beta decay can be viewed as the strong decay 
of neutron to n~p + system followed by the decay of the pion, we see that Eq. ( \L31\) is true under the 
same assumption as Eq.([70|). Namely, the pion should be much lighter than any other state with 
unnatural parity and zero baryon charge. Also the expectation value of the quark scalar operator 
averaged over pion is 

777*2 

X* = (7r*|gg|7r*) = . (132) 

m 

The other relations depend on the additional model assumptions. Starting with the ratio m*(p)/m 
we find in the straightforward generalization of NJL model 

= 44 • (133) 
m «(0) V ; 



This relation is referred to in the paper as Nambu scaling. The QCD sum rules prompt a more 
complicated dependence, presented by Eq.(|9T|) with the function 

F{p) = z l -7-rz— (134) 

l + av(p)/p 
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where a ~ —0.2 ]2"I| , pE3"| , — see also Sec. 5. Another assumption, expressed by Brown-Rho scaling 
equation (|85|) predicts a slower decrease of m*(p). Note, that Eq.(|85D is based on existence of a single 
length scale, while there are two at least: p F l and Aq^ d . 

The experimental situation with A-isobar mass in nuclear matter is not quite clear at the moment. 
The result on the total photon-nucleus cross section indicates that the mass m* A does not decrease 
in the medium [|110| , while the nucleon mass m*(p) diminishes with p. On the other hand, the 



experimental data for total pion-nucleus cross sections are consistent with the mass m* A decreasing in 
the matter As to calculations, the description within the Skyrmion model |93| predicts that m* A 

decreases in nuclear matter and m* A — m* < m^—m. Assuming the Additive Quark Model prediction 
for the scalar field-baryon couplings g S NN = 3sAA we come to the equation m* A — m* = — m. The 
Brown-Rho scaling leads to still smaller shift m* A — m* = [m*(m& — m)]/m. 

Now we present the results of the self-consistent calculations of the condensate under various 
assumptions on the dependence ^(^(p)) — Eq. (|130| ). 

In Fig. 7 we show the results with BR1 scaling of the nucleon mass — Eq.(^) for different sets of 
FFST parameters. Following [ |107|| we try the values which were obtained at saturation density p = po 
- set" a" defined in 4.4.1. We assume that they do not change with density. We use also another set 
of parameters = jna = 7aa = 0.7 (set "b"), presented in [|75j. The dependence on the behaviour 
m A.(p) appears to be more pronounced for set "a" of the FFST parameters. The calculations were 
carried out under the assumption /* = f n . Thus, the pion mass drops somewhat faster, than in BR1 
scaling with decreasing /*, in order to save in- medium GMOR relation — Eq. (|70|) . The nucleon 
effective mass at saturation density appears to be quenched somewhat less, than in QHD models, 
being closer to the value, preferred by FFST approaches |7I| |107 . 



Assuming that the pion mass does not change in medium we find strong dependence on the values 
of FFST parameters. For the choice "a" the self-consistent solution disappears before the density 



reaches the saturation value- see dotted curve "2" in Fig. 7. We explained in our paper ||103|| , how 
it happens technically. 

The results, obtained with Nambu scaling for the nucleon mass being assumed — Eq.( |133| ), are 
shown in Fig. 8. We put m* = and thus f* ~ Kp)| 1/2 , following GMOR. The three curves 
illustrate the dependence of the results on the assumption of the in-medium behaviour of the isobar 
mass. 

One of the results of this subsection is that the improved (self-consistent) approach excludes 
the possibility of the pion condensation at relatively small densities. On the other hand Dickhoff 
et al. | |101| | carried out self-consistent description of the particle-hole interactions by inclusion of 
the induced interactions to all orders. This shifts the point of the pion condensation to the higher 
densities. The rigorous analysis should include both aspects of self-consistency. 



4.4.4 Accumulation of isobars as a possible first phase transition 

While the density increases, the Fermi momentum and the energy of the nucleon at the Fermi surface 
increase too. At some value of p it becomes energetically favourable to start the formation of the 
Fermi sea of the baryons of another sort instead of adding new nucleons. This phase transition takes 
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Figure 7: The quark scalar condensate n(p) and nucleon effective mass m*(p), calculated in framework 
of BR1 scaling of the nucleon mass. The dash-double-dotted line shows the gas approximation. The 
dotted lines 1 and 2 present the results under assumption m*(p) = m % for the sets "b" and "a" of 
FFST parameters. The other curves present the results obtained under assumption f*(p) = f w and 
illustrate dependence on the choice of the values of FFST parameters and on the assumed behaviour 
m A (p). The solid and dashed curves are obtained for the set "a" of FFST parameters with the BR 
assumption m* A — m A = (m* — m) 1 ^ and for m* A = m A correspondingly. The two other curves are 
obtained for the set "b" under BR scaling assumption for the isobar mass (dot- dashed curve) and 
under assumption that the isobar mass does not change in medium (long-dashed curve). 
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Figure 8: The dependence n(p) and m*(p) under the assumption of Nambu scaling of the nucleon 
mass and m^p) = m n . The three curves illustrate the dependence on the behaviour of m* A (p). Solid 
line corresponds to BR scaling m A — = (m* — m) — , dashed curve to m* A — = m* — m while 
the dashed-dotted curve is obtained for m* A = m^. The calculations were carried out with the choice 
"b" of FFST parameters. Long-dashed line in Fig. 8a shows the gas approximation. 
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place at the value p a , determined by the condition 



(Pa) = (p 2 F + m* 2 (Pa)) 1/2 (135) 



rn 



B 



with pp a being the value of Fermi momentum, corresponding to p a , while m* B is the mass of the 
second lightest baryon at p = p a . 

The vacuum values of A and S + hyperon masses are respectively 115 MeV and 43 MeV smaller 
than that of A-isobar. However, both experimental and theoretical data confirm that the hyperons 
interact with the scalar fields much weaker than the nucleons. Thus, at least in framework of certain 
assumptions on the behaviour of m* A (p), the hierarchy of the baryon masses changes in medium. 
(The investigations of the problem are devoted mostly to the case of neutron or strongly asymmetric 
matter because of the astrophysical applications. See, however, the paper of Pandharipande ||112|| ). 



The delta isobar can become the second lightest baryon state. In this case accumulation of A-isobars 
in the ground state is the first phase transition in nuclear matter. Such possibility was considered in 
several papers ||113|| - ||il~6l1 . 



Under the assumption of BR1 scaling the accumulation of A-isobar takes place at p a ~ 3po, being 



the first phase transition. The value of p a is consistent with the result of Boguta [114 . 



5 QCD sum rules 

5.1 QCD sum rules in vacuum 

Here we review briefly the main ideas of the method. There are several detailed reviews on the 
subject — see, e.g., [ |117|| . Here we focus on the points, which will be needed for the application of 



the approach to the case of nuclear matter. 

The main idea is to establish a correspondence between descriptions of the function G, introduced 
in Subsec.3.6 in terms of hadronic and quark-gluon degrees of freedom. (Recall that G describes the 
propagation of the system with the quantum numbers of the nucleon). The method is based on the 
fundamental feature of QCD, known as the asymptotic freedom. This means, that at q 2 — > — oo the 
function G(q 2 ) can be presented as the power series of q~ 2 and QCD coupling a s . The coefficients 
of the expansion are the expectation values of local operators constructed of quark and gluon fields, 
which are called "condensates" . Thus such presentation, known as operator product expansion (OPE) 
|118| |, provides the perturbative expansion of short-distance effects, while all the nonperturbative 



physics is contained in the condensates. 

The correspondence between the hadron and quark-gluon descriptions is based on Eq.(^). The 
empirical data are used for the spectral function ImG(A; 2 ) in the rhs of Eq. fl93|) . Namely, we know, 
that the lowest lying state is the bound state of three quarks, which manifests itself as a pole in the 
(unknown) point k 2 = m 2 . Assuming, that the next singularity is the branching point k 2 = W 2 h = 
(m + m n ) 2 , one can write exact presentation 



Im G(k 2 ) = \ 2 5{k 2 - m 2 ) + f{k 2 )0{k 2 -W 2 h ) (136) 



with A 2 being the residue at the pole. Substituting rhs of Eq. ( |136|) into Eq. (^) and employing q 



v-2 
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power expansion in lhs, i.e. putting 



G{q 2 ) = G OPE {q 2 ) (137) 
one finds certain connections between quark-gluon and hadron presentations 

G P E (q 2 ) = + -] dk " ■ ( 138 ) 

m z — q z 7r J k z — q z 

pn, 

Of course, the detailed structure of the spectral density f{k 2 ) cannot be resolved in such approach. 
The further approximations can be prompted by asymptotic behaviour 

f{k 2 ) = ±AG PE(k 2 ) (139) 

at k 2 ^> \q 2 \ with A denoting the discontinuity. The discontinuity is caused by the logarithmic 
contributions of the perturbative OPE terms. The usual ansatz consist in extrapolation of Eq.( |139| ) 
to the lower values of k 2 , replacing also the physical threshold W 2 h by the unknown effective threshold 
W 2 } i.e. 



7r ■/ k 2 — q 2 2ni J ^ k 2 — q 

and thus 



w\ w 2 

pn 



G ° Mq } = ^q~ 2 + 2^ J k 2 -q 2 dk ■ (M1) 

w 2 

The lhs of Eq.( |141|) contains QCD condensates. The rhs of Eq. (|141[) contains three unknown pa- 
rameters: m, X 2 and W 2 . Of course, Eq. ( |141| ) makes sense only if the first term of the rhs, treated 
exactly is larger than the second term, treated approximately. 

The approximation G(q 2 ) ~ GoPE{q 2 ) becomes increasingly true while the value \q 2 \ increases. 
On the contrary, the "pole+continuum" model in the rhs of Eq. ( |141| ) becomes more accurate while 
|g 2 | decreases. The analytical dependence of the lhs and rhs of Eq. ([L41]) on q 2 is quite different. The 
important assumption is that they are close in certain intermediate region of the values of q 2 , being 
close also to the true function G(q 2 ). 

To improve the overlap of the QCD and phenomenological descriptions, one usually applies the 
Borel transform, defined as 



Bf(Q 2 ) = Jim (-JL ) /( g 2 ) = f(M 2 ) (142) 



Q 2 = -q 2 - M 2 = Q 2 /n 
with M called the Borel mass. There are several useful features of the Borel transform. 



1. It removes the divergent terms in the lhs of Eqs. ( |138| ) and (|141|) which are caused by the free 
quark loops. This happens, since the Borel transform eliminates all the polynomials in q 2 . 
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2. It emphasise the contribution of the lowest lying states in rhs of Eq. (|141 ) due to the relation 



B 



Q2 _|_ m 2 



-m 2 /M 2 



(143) 



3. It improves the OPE series, since 



B 



(Q) 



2\-n 



(n-1)! 



(M 



2\l-n 



(144) 



Applying Borel transform to both sides of Eq.( 141 ) one finds 

oo 

G OPE (M 2 ) = A 2 e- m2 / M2 + -i, / dk'e'^ 2 ■ AG OPE (k 2 ) . 

2m J 



(145) 



w 2 



Such relations are known as QCD sum rules. If both rhs and lhs of Eq. (|141[) were calculated exactly, 
the relation would be independent on M 2 . However, certain approximations are made in both sides. 
The basic assumption is that there exists a range of M 2 for which the two sides have a good overlap, 
approximating also the true function G(M 2 ). 

The lhs of Eq. ( |145 ) can be obtained by presenting the function G(q 2 ), which is often called 
"correlation function" or "correlator" as (strictly speaking, G depends on the components of vector 
q also through the trivial term q) 



G{q 2 ) = i [ d 4 xe i{qx) (0\T{ri(x)fj(0)}\0) 



(146) 



with 7] being the local operator with the proton quantum numbers. It was shown in ||1 19|| that there 
are three independent operators r\ 



Vi = 



(u T a C lll u b ) l5 ^d c -e abc J r, 2 = (u T a Ca, u u b )a^^d c 



.abc 



[ulCj^u^dc - {u T a C^^d b )^u c 



.abc 



(147) 



where T denotes the transpose in Dirac space and C is the charge conjugation matrix. However, the 



operator r] 2 provides strong admixture of the states with negative parity [|119|| . As to the operator 773, 
it provides large contribution of the states with spin 3/2 |119| . Thus, the calculations with f] = rji 
are most convincing. We shall assume rj = r]\ in the further analysis. 



The correlation function has the form 



G{q) = G q {q 2 )-q + G s {q 2 )-I 



(148) 



with / standing for the unit 4x4 matrix. The leading OPE contribution to G q comes from the 
loop with three free quarks. If the quark masses m u>d are neglected, the leading OPE term in G s 
comes from the exchange by the quarks between the system described by operator 77 and vacuum. 
Technically this means, that the contribution comes from the second term of the quark propagator 
in vacuum 



(01^(^(0)10) 



X a /3 



1 



2tt 2 x 4 4 



Er^(o|gr A g|o) + o( 



x 



(149) 
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Figure 9: The Feynman diagrams, describing the lowest order OPE contribution to the nucleon 
correlator in vacuum. The helix line stands for the system with the quantum numbers of proton. 
The solid lines denote the quarks. The light circles denote the vacuum expectation value. 

where only the contribution with A = 1 (see Eq.(|l^)) survives. This is illustrated by Fig.9. The 
higher order terms come from exchange by soft gluons between vacuum and free quarks carrying 
hard momenta. Next comes the four-quark condensate which can be viewed as the expansion of the 
two-quark propagator, similar to Eq.( |149| ). 

Direct calculation provides for massless quarks [|98[1 



qOPE 
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h 



(150) 
(151) 



with the condensates g(0) = (0|— G 2 |0), k(0) = (0|gg|0) and h = (0\uuuu\0). The terms, containing 
polynomials of q 2 are omitted, since they will be eliminated by the Borel transform. This leads to 



the sum rules d 



„ W 2 \ 1 9 W 2 ' 
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1 Co = \ 2 e- m2 ' M2 
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y2 e -m 2 /M 2 



(152) 
(153) 



with traditional notations a = — 27t 2 k(0), b = (27r) 2 g(0), A 2 = 327r 4 A 2 



Ea(x) 



1 - e- x 



EAx) 



1 - (l + x)e 



Eo(x) 



1 - 



X 



+ x+l]e~ 



Also Co = (27r) ho- Here we omitted the anomalous dimensions, which account for the most im- 
portant corrections of the order a s , enhanced by the "large logarithms". The radiative corrections 



were shown to provide smaller contributions, as well as the higher order power corrections p8fl . The 
matching of the lhs and rhs of Eqs.(152),(153) was found in |)8| for the domain 



0.8 GeV 2 < M < 1.4 GeV 2 



(154) 
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As one can see from Eg. Ql53| ), the nucleon mass turns to zero if (0|gg|0) = 0. Hence, the mass is 
determined by the exchange by quarks between the our system and vacuum. 

The method was applied successfully to calculation of the static characteristics of the nucleons 
reproducing the values of its mass [[39], |97|, [98|] as well as of magnetic moment @] and of the axial 
coupling constant 
approach ||122[| 



121 1 . The proton structure functions also were analysed in framework of the 



5.2 Proton dynamics in nuclear matter 

Now we extent the sum-rule approach to the investigation of the characteristics of the proton in 
nuclear matter. The extension is not straightforward. This is mostly because the spectrum of 
correlation function in medium 

G m {q) = i J d 4 xe^ qx \M\T{ri(x)fj(0)}\M) (155) 

is much more complicated, than that of the vacuum correlator G(q 2 ). The singularities of the 
correlator can be connected with the proton placed into the matter, as well as with the matter itself. 
One of the problems is to find the proper variables, which would enable us to focus on the properties 
of our probe proton. 



5.2.1 Choice of the variables 

Searching for the analogy in the earlier investigations one can find two different approaches. Basing 
on the analogy with the QCD sum rules in vacuum, one should build the dispersion relation in the 
variable q 2 . The physical meaning of the shift of the position of the proton pole is expressed by 



Eq.(^). Another analogy is the Lehmann representation ||120| , which is dispersion relation for the 
nucleon propagator in medium gN(qo, i n the time component go- Such dispersion relation would 
contain all possible excited states of the matter in rhs. Thus, we expect the dispersion relations in 
q 2 to be a more reasonable choice in our case. 

It is instructive to adduce the propagation of the photon with the energy u and three dimensional 
momentum k in medium. The vacuum propagator is D 1 ~ [uj 2 — k 2 ]~ l . Being considered as the 
function of q 2 = uo 2 — k 2 it has a pole at q 2 = 0. The propagator in medium is D™ ~ [u 2 e(uj) — k 2 }^ 1 . 
The dielectric function e(u) depends on the structure of the matter, making .D™(u;) a complicated 
function. However, the function D™(q 2 ) still has a simple pole, shifted to the value q 2 ^ = uj 2 {1 — s{oj)). 
A straightforward calculation of the new value q^ is a complicated problem. The same refers to the 
proton in-medium. The sum rules are expected to provide the value in some indirect way. 

Thus, we try to build the dispersion relations in q 2 . Since the Lorentz invariance is lost, the 
correlator G m (q) depends on two variables. Considering the matter as the system of A nucleons with 
momenta Pi, introduce vector 

P = ^ , (156) 

which is thus p ~ (m, 0) in the rest frame of the matter. The correlator can be presented as 
G m (q) = G m (q 2 ,(p(p,q)) with the arbitrary choice of the function (p(p,q), which is kept constant in 
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the dispersion relations. This is rather formal statement, and there should be physical reasons for 
the choice. 

To make the proper choice of the function (p(p, q), let us consider the matrix element, which enters 
Eq.flliD 

(M\T{ v (x)f)(0)}\M) = (M A \ V (x)\M A+1 )(M A+1 \fj(0)\M A )e{x Q ) - 

- (M A \fj(0)\M A _ l )(M A _. 1 \ r} (x)\M A ) 6(-x ) (157) 

with | Ma) standing for the ground state of the matter, while |Ma±i) are the systems with baryon 
numbers A ± 1. The summation over these states is implied. The matrix element (Ma+i|?7|Ma) 
contains the term (N\r]\0) which adds the nucleon to the Fermi surface of the state |Ma). If the 
interactions between this nucleon and the other ones are neglected, it is just the pole at q 2 = m 2 . 
Now we include the interactions. The amplitudes of the nucleon interactions with the nucleons of 
the matter are known to have singularities in variables Si = (pi + q) 2 . These singularities correspond 
to excitation of two nucleons in the state |Ma+i). Thus, they are connected with the properties of 
the matter itself. To avoid these singularities we fix 

<p(p,q) = s = maxsj = 4i? 2 F (158) 

with Eqf being the relativistic value of nucleon energy at the Fermi surface. Neglecting the terms of 
the order p 2 F /m 2 we can assume pi = (m, 0) and thus 

s = 4m 2 . (159) 

Our choice of the value of s corresponds to |g| = pf (in the simplified case, expressed by Eq.( |159| ) 
\q\ =0). However, varying the value of s we can find the position of the nucleon poles, corresponding 
to other values of \q\. 

Let us look at what happens to the nucleon pole q 2 = m 2 after we included the interactions with 
the matter. The self-energy insertions £ modify the free nucleon propagator g% to gN with 

(9n)- 1 = (g° N r l -z (160) 

— see Fig. 10. 

In the mean field approximation (Fig. 10a) the function E does not contain additional intermediate 
states. It does not cause additional singularities in the correlator G m (q 2 , s). The position of the pole 
is just shifted by the value, which does not depend on s. (Note that this does not mean that in the 
mean field approximation the condition s = const can be dropped. Some other contributions to the 
matrix element (Ma+i|^|Ma) are singular in s. Say, there is the term (B\fj\0) with B standing for 
the system, containing the nucleon and mesons. If the mesons are absorbed by the state \M A ), we 
come to the box diagram (Fig. 11) with the branching point, starting at s — 4m 2 ). 

Leaving the framework of the mean-field approximation we find Hartree self-energy diagrams 
(Fig. 10b) depending on s. The latter is kept constant in our approach. Hence, no additional singu- 
larities emerge in this case as well. 

The situation becomes more complicated if we take into account the Fock (exchange) diagrams 
(Fig. 10c). The self-energy insertions depend on the variable u — (p — q) 2 . The contribution of these 
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Figure 10: Self-energy insertions to the nucleon pole contribution to the correlator. The solid 
lines denote the nucleons, the helix line stands for the correlator. Fig. 10a corresponds to mean- field 
approximation. Fig. 10b shows the self-energy in the direct channel. Fig. 10c presents the exchange 
contribution. 




Figure 11: One of the contributions, providing the branching point at s — Am 2 . The bold solid line 
denotes the nucleon of the matter. The dashed lines stands for the meson systems. 
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terms shifts the nucleon pole and gives birth to additional singularities corresponding to real states 
with baryon number equal to zero. They are the poles at the points u = m 2 , with m x denoting the 
masses of the mesons (tt, lo, etc.), and the cuts running to the right from the point q 2 = m 2 + 2m 2 . 
The latter value is the position of the branching point corresponding to the real two-pion state in 
the w-channel. 

Thus the single-nucleon states |-B±i) cause the pole q 2 = m^, a set of poles corresponding to 
the states with baryon number B = and a set of branching points. The lowest-lying one is 
q 2 = m 2 + 2m 2 . Note that the antinucleon corresponding to go = ~ m generates the pole q 2 = 5m 2 
shifted far to the right from the lowest-lying one. 

The lowest-lying branching point q 2 = m 2 +2m 2 is separated from the position of the pole q 2 = m 2 
by a much smaller distance than in the case of the vacuum (q 2 = m 2 + 2mm 7T in the latter case). Note, 
however, that at the very threshold the discounting is quenched since the vertices contain moments 
of the intermediate pions. Thus the branching points can be considered as a separated from the pole 
q 2 = m 2 . Note also that for the same reason the residue at the pole q 2 = m 2 + 2m 2 in the ^-channel 
vanishes. 

As it was shown in the work fll3|| , all the other singularities of the correlator G m (q 2 , s) in q 2 are 
lying to the right from the nucleon pole until we include the three-nucleon terms. Thus, they are 
accounted for by the continuum and suppressed by Borel transform. To prove the dispersion relation 
we must be sure of the possibility of contour integration in the complex q 2 plane. This cannot be 
done on an axiomatic level. However a strong argument in support of the possibility is the analytical 
continuation from the region of large real q 2 . At these values of q 2 the asymptotic freedom of QCD 
enables one to find an explicit expression of the integrand. The integral over large circle gives a 
non-vanishing contribution. However, the latter contains only the finite polynomials in q 2 which are 
eliminated by the Borel transform. 

Thus we expect "pole+continuum" model to be valid for the spectrum of the correlator G m (q 2 , s). 

The situation becomes more complicated if we include the three-nucleon interactions 155). The 



probe proton, created by the operator rj can interact with n nucleons of the matter. The corresponding 
amplitudes depend on the variable s n = (np + q) 2 . For n>2 this causes the cuts, running to the left 
from the point q 2 = m 2 . This requires somewhat more complicated model of the spectrum. From the 
point of view of expansion in powers of p this means, that "pole + continuum" model is legitimate 
until the terms of the order p 2 are included. 

5.2.2 Operator expansion 

Following our general strategy, we shall try to obtain the leading terms of expansion of the correlator 

G m (q) = G™(q)q + G™(q)p + G™{q)I , (161) 

in powers of q~ 2 . Note, that the condition s = const, which we needed for separation of the singu- 
larities, connected with our probe proton, provides 

-> const (162) 
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at q 2 — > — oo. This is just the condition which insures the operator expansion in deep- inelastic 
scattering (see, e.g., the book of Ioffe et al. ||123|| ). It is not necessary in our case. However, 



the physical meaning of some of the condensates, say, that of ((p a (a)a n ) — Eq . (|50|) becomes most 
transparent in this very kinematics. 

The problem is more complicated than in vacuum, since each of the terms of the expansion in 
powers of q 2 provides, generally speaking, infinite number of the condensates. Present each of the 
components G™ (i = q,p, s) of the correlator G m 

Gf = [ d^xS^Ti{x) . (163) 



The function Tj(x) contains in- medium expectation values of the products of QCD operators in 
space-time points "0" and "x" with an operator at the point x defined by Eq.(fl^). Each in-medium 
expectation value, containing covariant derivative is proportional to the vector p^. This can be 
easily generalized for the case of the larger number of derivatives. Thus the correlators take the form 
G™ — J2 n C n {pV q ) n fi(q 2 ). For the contributions fi{q 2 ) ~ {q 2 )~ k the terms {pV q ) n fi(q 2 ) are of the 
same order. This is the "price" for the choice of kinematics s = const. Fortunately, the leading terms 
of the operator expansion contain the logarithmic loops and thus can be expressed through the finite 
number of the condensates Q . 

The leading terms of the operator expansion can be obtained by replacing the free quark propa- 
gators by those in medium 

(M\TM*)M0)M = ^|-E^ r i( M l^(0)r A ^(x)|M) (164) 

with the matrices T A being defined by Eq. (JIB|) . Operator ip(x) is defined by Eq.fl4"5]). While looking 
for the lowest order term of the operator expansion we can put x 2 = in the second term of the rhs 
of Eq.( 164 ). In the sum over A the contributions with A = 3, 4 vanish due to the parity conservation 



by strong interactions, the one with A = 5 turns to zero in any uniform system. Thus, only the terms 
with A — 1, 2 survive. Looking for the lowest order density effects, we assume that propagation of 
one of the quarks of the correlator G m is influenced by the medium. Hence, the term with A = 1 
contributes to the scalar structure G™, while that with A = 2 — to the vector structures G™ and 

Op . 

G™ = ^q 2 \n(-q 2 )K(p) (165) 

G? = (1 2 ) 2 M-<1 2 ) + -^ (s-m 2 -q 2 )\n(-q 2 )v(p) (166) 

G™ = ^q 2 \n(-q 2 )v(p). (167) 

Thus the correlator G™ can be just obtained from the vacuum correlator G s by replacing of k;(0) 
by n(p). The correlator G™ obtains additional contribution proportional to the vector condensate 
v(p). Also, the correlator G™, which vanishes in vacuum is proportional to v(p). These terms are 
illustrated by Fig. 12 a,b. 

Turn now to the next OPE terms. Start with the structure G™. In the case of vacuum there is a 
contribution which behaves as ln(— q 2 ), which is proportional to the condensate (Oltp^G^a^ ^ip\0) . 
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Figure 12: The Feynman diagrams, contributing to the leading terms of OPE of the nucleon correlator 
in medium. The dark blob denotes in-medium expectation values. The dotted lines stand for gluons. 
The others notations coincide with those of Fig.9. 



However, similar term comes from expansion of expectation value (0|g(0)g(x)|0) in powers of x 2 . 
The two terms cancel |39|. Similar cancellation takes place in medium ||99|| . However there is a 
contribution, caused by the second term of rhs of Eq.([l5|). It does not vanish identically, but it can 
be neglected due to Eq. ([54"D . Hence, the next OPE term in rhs of Eq.( |165| ) can be obtained by simple 



replacement of the condensates k(0) and g(0) in the second term of Eq. (|151|) by n(p) and g(p)- 

The next-to-leading order corrections to the correlators G™ p come from expansion of the expec- 
tation value (M|-0(O)7oV'(^)|^)- In the lowest order of x 2 expansion the matrix element can be 
presented through the moments of the deep-inelastic scattering (DIS) nucleon structure functions - 
Eqs. ( |4"9| ) and (|50|). Since the medium effects in DIS are known to be small we limit ourselves to the 
gas approximation at this point. The main contributions to q~ 2 expansion compose the series of the 
terms [(s — m 2 ) / q 2 ] n (a n ) with (a n ) denoting n-th moment of the structure function. Being expressed 



in a closed form, they change Eqs. ( |166|) and ( |167| ) to 
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(168) 
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J da(f)b(a) ln(g — pa) 2 ■ p 
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± 

°p = ^ q2 J daF ^ - pa ^ 2 ■ p ( 169 ) 



with F(a) being structure function, normalized as Jq 1 daF(a) = 3. The function F(a) can be 
presented also as F(a) = 4> a (a) with (fi a defined by Eq.flSop. Another leading OPE term is caused 
by modification of the value of gluon condensate. This is expressed by changing of the value g(0) in 
the second term in rhs of Eq.( 150 ) to g(p) - Fig. 12b. 

The higher order OPE terms lead to the contributions which decrease as q~ 2 . One of them is 
caused by the lowest order power correction to the first moment of the structure function. This 
term is expressed through the factor £ which is determined by Eq.(51), being calculated in ||42|| . The 
other corrections of this order come from the four-quark condensates Q AB , defined by Eq . fl57|) . The 
correlator G™ contains the condensate Q 12 . The vector part includes the condensates Q 11 and Q 22 . 
Another contribution of this order comes from the replacement of the condensates g(0) and k(0) in 
the last term in rhs of Eq.(151) by their in-medium values - Fig. 12c, d. 



5.2.3 Building up the sum rules 

To construct the rhs of the sum rules, consider the nucleon propagator = (H — E)^ 1 with 
E standing for the nucleon energy, while the Hamiltonian H in the mean field approximation is 
presented by Eq . (|73|) . Beyond the mean field approximation the potentials and $ should be 
replaced by vector and scalar self-energies 

V, = El; El=p^ + q^; $ = . (170) 

Thus, under condition s = 4m 2 — Eq.( |159| ) 



g(l-S*)-pSP + m + S* 

9n = Z 171 

q 2 -m 2 m 

with 

m m = m + U ; U = m(S 9 + S p ) + S s , (172) 

while 



1 



(173) 



(1 - E«)(l - + Ep) ' 
Of course, S 9 = in mean field approximation. 

The Borel-transformed sum rules for in-medium correlators in the assumed "pole+continuum" 
model for the spectrum are: 

£™(M 2 ) = \ 2 m e- m ™l M \l-Y?) (174) 
£-(M 2 ) = -\ 2 m e- m ™ /M2 Y? (175) 
£™(M 2 ) = \ 2 m e- m ™ /M \m + i: s ) (176) 
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with A 2 ^ = 32n 4: X 2 n Z with A^ standing for the residue in the nucleon pole (see similar definition in 
vacuum — Eqs.(p2D and (pSj)). The lhs of Eqs. (|l74l)-(|17^) are: 



£™(M 2 ) = M% (^) L~*/» - ^ 



x (F/x(a)) - 2m 2 M 2 E 



'W 2 \ fw 2 \ 

• (Fp(a)a) + Am 2 M 2 E (<f> bf z(a)) 



M 2 



M 2 
4 



1 \ M 2 ) 



(W 2 \ 3 

+n 2 M 2 E s(p) + i m 2 ( S - m 2 )(# a/ i> + 2m 4 (M + 3 (2vr) 4 g^(p)L 



£™(M 2 ^ 



3 2 /TI/2 \ 



M 2 



(F/x(a)) P L- 



2\2 



£™(M 2 ) = (2vr 2 ) 



M 4 E ( W| 
1 1 M 2 



(27T) 

12 



(2vr) 4 Q^ 



(177) 
(178) 
(179) 



In Eqs. ( |17?1 )- (|179|) (-0) = Jq daijj(a) for any function ip, F is the structure function, the function 

p{a) = e W ( - {s -^ )a + m2a2 ) (180) 



*\ M 2 (l + a 
takes into account the terms [(s — m 2 )/M 2 ] n . The factor 

lnM 2 /A 2 



In v 2 /k 2 



181) 



accounts for the anomalous dimensions. Here A = 0.15 GeV is the QCD parameter while v = 0.5 GeV 
is the normalization point of the characteristics involved. 

Recall, that "pole+continuum" model is true until we do not touch the terms of the order p 2 . 
Thus, in the sum rules for the difference between in-medium and vacuum correlators we must limit 
ourselves to linear shifts of the parameters 



AC^M 2 ) = AV-VM* ^) - ^ exp ( AW 2 



A 2 



M 2 J 2L 4 / 9 



M 2 



£™(M 2 



^2^-m 2 /M 2 YP 



A£™(M 2 ) = A 2 e" m2 / M2 (m ^ + £ s - ^^-) - 2aW exp 



A 2 



M 2 



M 2 



AW Z 



(182) 
(183) 
(184) 



Here A denotes the difference between in-medium and vacuum values. However the self-energy T, q 
and S s cannot be determined separately, since only the sum + S s can be extracted. 

m can be obtained: using Eq.( |172j ) 



Anyway, the shift of the position of the nucleon pole m 
we find 



A£™ - m£™ - mA£; 



AmA e 



2 -m 2 /M 2 



+ W z e 



2 -m 2 /M 2 



W 
l2L 4 / 9 



m - 2a) AW 1 



185) 
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bmce the value (W 2 /2L 4 / 9 )m - 2a is numerically small, one can write approximate sum rule, ne- 
glecting the second term in rhs of Eq. ( |185| ) 

m 2 /M 2 

U = — — (A£™ - m£™ - mA£") , (186) 
or, assuming the sum rules in vacuum to be perfect 

e m 2 /M 2 

U = — ^— (C s - mC p - mC q ) (187) 
with the vacuum part cancelling exactly. 

2 

The two lowest order OPE terms ( without perturbative expansion in parameter ^rpr ) are 
presented by the first two terms of rhs of Eqs. ( |177|) , ( 179 ) and by the first term of rhs of Eq. ( |178|) . 



They are expressed through the condensates v(p), k(p) and g(p) and through the moments of the 
nucleon functions 4> a ^ introduced in Subsec. 2.6. The values of the lowest moments of the structure 
function F(a) = <fia{ot) are well known from experimental data. By using the value of £ a and 
employing relations, presented by Eq.([53|) one can find the lowest moments of the function Only 
the first moment of the function 0^ and thus the first and second moments of the function cf) a appeared 
to be numerically important. Thus, at least in the gas approximation all the contributions to the lhs 
of the sum rules can be either calculated in the model-independent way or related to the observables 



43 . The scalar condensate is the most important parameter beyond the gas approximation pO 



The model calculations have been carried out in this case. 

The next order of OPE includes explicitly the moments of the functions 9 a ^ defined in Subsec. 2. 6. 
It includes also the four-quark condensates Q n , Q 12 and Q 22 . Using Eqs. (|53|) one can find that only 
the first moment of the function 9 a is numerically important while the moments of the function 0& 
can be neglected. The condensate Q 12 can be obtained easily by using Eq . (|58|) . The uncertainties 
of the values of the other four-quark condensates Q u is the main obstacle for decisive quantitative 
predictions, based on Eqs. ( |182j )- (|187|) . The scalar four-quark condensate Q 11 may appear to be a 
challenge for the convergence of OPE due to the large value of the second term in rhs of Eq.(^). This 
may be a signal that large numbers are involved. Fortunately, the only calculation of Q 11 carried out 



in [46| demonstrated that there is a large cancellation between the model-dependent first term in rhs 
of Eq.fl5"9|) and the second one, which is to large extent model-independent. However, assuming the 
result presented in we still find this contribution to be numerically important. 

We can try (at least for illustrative reasons) to get rid of this term in two ways. One of them is 
to ignore its contribution. The reason is that it corresponds to exchange by a quark system with the 
quantum numbers of a scalar channel between our probe proton and the matter. On the other hand, 
it contributes to the vector structure of the correlator G™, and thus to the vector structure q of the 
propagator of the nucleon with the momentum q. Such terms are not forbidden by any physical law. 
However most of QHD calculations are successful without such contributions. Thus the appearance 
of the terms with such structure, having a noticeable magnitude is unlikely. (Of course, this is not a 
physical argument, but rather an excuse for trying this version). The other possibility is to eliminate 
the contribution by calculation of the derivative with respect to M 2 . The two ways provide relatively 
close results. 
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5.2.4 The structure of the potential energy 



Under the conditions, described above, we find that the rhs of Eq.( |186|) is a slowly varying function 
of M 2 in the interval, defined by Eq. ( |154| ). Among the moments of the structure function the two 
first ones appeared to be numerically important. Thus we find 



U( P ) 



66v(p) + 70v 2 (p) - 



10Ag(p) 



- 32 Ak(p) 



GeV 



-2 



Here v(p) = 3p is the vector condensate — Eq. (|T6|) , An(p) = /t(p) — k(0) is the in-medium change 
of the scalar condensate — Eq.(62). The condensate t>2(p), determined as 



(M$^D„iJ)\M) = [g^ 



p 2 



0.3p. 



is connected to the second moment of the nucleon structure function. Numerically V2(p) 
Finally, Ag(p) is the shift of the gluon condensate, expressed by Eq . (|39|) . 

Thus the problem of presenting the nucleon potential energy through in-medium condensates is 
solved. At the saturation value p = po we find U = —36 MeV in the gas approximation. This should 
be considered as a satisfactory result for such a rough model. This is increasingly true, since there 
is a compensation of large positive and negative values in rhs of Eq.( |186| ). 

Note that the simplest account of nonlinear terms signals on the possible saturation mechanism. 
Following the discussion of Subsection 2.7 and assuming the chiral limit, present Ak(p) = S/m — 
3.2(p F /p F0 )p. Thus we obtain the potential 



U(p) 



198 -42—) 

mj po 



133 



4/3' 



MeV 



(189) 



After adding the kinetic energy it provides the minimum of the functional S(p) defined by Eq. (f?9|) 
at S = 62.8 MeV, which is consistent with experimental data — Eq. fl2"5|) . The binding energy 
is £ = —9 MeV. The incompressibility coefficient K = 9po(d 2 e/dp 2 ) also has a reasonable value 
K = 182 MeV. 

Of course, the results for the saturation should not be taken too seriously. As we have seen in 
Sect. 3, the structure of the nonlinear terms of the condensate is much more complicated. Also, the 
result is very sensitive to the exact value of S-term. Say, assuming it to be larger by the magnitude 
of 2 MeV, we find the Fermi momentum at the saturation point about 1/3 larger than p F0 . Thus, 
the value of the saturation density becomes about 2.5 times larger than p . Such sharp dependence 
is caused by the form of the nonlinear term in the potential energy equation — Eq.( |189| ). The form 
of the term is due to oversimplified treatment of nonlinear effects. However the result can be the 
sign, that further development of the approach may appear to be fruitful. 



5.2.5 Relation to conventional models and new points 

We obtained a simple mechanism of formation of the potential energy. Recall that Ioffe analysis of 
QCD vacuum sum rules |)7|| provided the mechanism of formation of nucleon mass as due to the 
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Figure 13: Contribution to the scalar structure of the correlator from the exchange of qq pairs with 
quantum numbers of vector mesons between the correlator and the matter. The dark blob stands for 
the vector condensate. The light circles denote the vacuum expectation value. The generated term 
is unusual for QHD. 



exchange by quarks between the probe nucleon and the quark-antiquark pairs of vacuum. In the 
nuclear matter the new mass is formed by the exchange with the modified distribution of the quark- 
antiquark pairs and with the valence quarks. The modified distribution of qq pairs is described by 
the condensate n(p). At p close to p the modification is mostly due to the difference of the densities 
of qq pairs inside the free nucleons and in the free space. Similar exchange with the valence quarks 
is determined by the vector condensate v(p) and is described by the first term of rhs of Eq. (|188|) . 
The second term describes additional interaction which takes place during such exchange. These 
exchanges cause the shift of the position of the pole m m — m. While the interactions of the nucleons 
depend on the condensates An(p) and v(p), these condensates emerge due to the presence of the 
nucleons. Also, the nonlinear part of An is determined by NN interactions. Thus, there is certain 
analogy between QCD sum rules picture and NJL mechanism. 

As we have seen, the QCD sum rules can be viewed as connection between exchange by uncorre- 
cted qq pairs and exchange by strongly correlated pairs with the same quantum numbers (mesons). 
This results in connection between the Lorentz structures of correlators and in-medium nucleon prop- 
agators. In the leading terms of OPE the vector (scalar) structure is determined by vector (scalar) 
condensate. The large values (of about 250-300 MeV) of the first and the fourth terms in rhs of 
Eq. ( |188| ) provide thus the direct analogy with QHD picture. 

Note, however, that the sum rules, presented by Eqs. (|177|) - (|179|) contain also the terms, which 



are unusual for QHD approach. Indeed, the term Q 11 in Eq.( |177| ) enters the vector structure of the 
correlator (and thus, of the propagator of the nucleon) corresponding, however, to exchange by the 
vacuum quantum numbers with the matter. On the other hand, the last term of rhs of Eq . ( |1 79| ) 
treated in the gas approximation, corresponds to exchange by the quantum numbers of vector mesons. 
However, it appears in the scalar structure. This term originated from the four-quark condensate Q 12 . 
While the exact value of the condensate Q n is still obscure, the condensate Q 12 is easily calculated. 
This OPE term is shown in Fig. 13. It provides a noticeable contribution. 

Such terms do not emerge in the mean-filed approximation of QHD. They can be originated by 
more complicated structure of the nucleon-meson vertices. (Note that if the nucleons interacted 
through the four-fermion interaction, such terms would have emerged from the exchange interaction 
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124 . In most 



due to Fierz transform). 

Another approach, developed by Maryland group, was reviewed by Cohen et al. 
of the papers (except ||125| |) the Lehmann representation was a departure point. In framework of 
this approach the authors investigated the Lorentz structure of QCD sum rules |126| . They analysed 
detaily the dependence of the self-energies on the in-medium value of the scalar four-quark condensate 
|127|| . The approach was used for investigation of hyperons in nuclear matter ||128| . 

The approach is based on the dispersion relations in the time component go at fixed three- 
dimensional momentum |g|. It is not clear, if in this case the singularities, connected with the probe 
proton are separated from those of the matter itself. The fixed value of three-dimensional momenta 
is a proper characteristics for in-medium nucleon. That is why this was the choice of variables in 
Lehmann dispersion relation with the Fermi energy as a typical scale. The sum rules are dispersion 
relations rather for the correlation function with the possible states N, iV+pions, N*, etc. The scale 
of the energy is a different one and it is not clear, if this choice of variables is reasonable for QCD 
sum rules. 



5.3 Charge -symmetry breaking phenomena 
5.3.1 Nolen-Schiffer anomaly 

The nuclei consisting of equal numbers of protons and neutrons with one more proton or neutron 
added are known as the mirror nuclei. If the charge symmetry (known also as isospin symmetry) 
of strong interactions is assumed, the binding energy difference of mirror nuclei is determined by 
electromagnetic interactions only, the main contribution being caused by the interaction of the odd 
nucleon. Nolen and Schiffer ||129|| found the discrepancy between the experimental data and theoret- 



ical results on the electromagnetic contribution to the energy difference. This discrepancy appeared 
to be a growing function of atomic number A. It reaches the value of about 0.5 MeV at A = 40. 
Later the effect became known as Nolen-Schiffer anomaly (NSA). 

The NSA stimulated more detailed analysis of electromagnetic interactions in such systems. Auer- 
bach et al. |130|| studied the influence of Coulomb forces on core polarization. However, this did not 



explain the NSA. Bulgac and Shaginyan ||131|| attributed the whole NSA phenomena to the influ- 
ence of the nuclear surface on the electromagnetic interactions. Thus they predict NSA to vanish in 
infinite medium. 

However most of the publications on the subject contain the attempts to explain NSA by the 
charge symmetry breaking (CSB) by the strong interactions at the hadronic level. The CSB potentials 
of NN interactions were reviewed by Miller et al. ||132 |. Some of phenomenological potentials 



described the NSA, but contradicted the experimental data on CSB effects in NN scattering. The 
meson-exchange potentials contain CSB effects by inclusion of p — u mixing. This explains the large 
part of NSA, but not the whole effect. 

On the quark level the CSB effects in the strong interactions are due to the nonzero value of the 
difference of the quark masses 

m d « 7 MeV ; m u « 4 MeV ; \i = m d - m u « 3 MeV . (190) 
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Several quark models have been used for investigation of NSA by calculation of neutron-proton mass 
difference in nuclear matter (recall that in vacuum m n — m p = 1.8 MeV, while the Coulomb energy 
difference is -0.5 MeV. Hence, the shift caused by strong interactions is 5m = 2.3 MeV). Henley and 



Krein ||133|| used the NJL model for the quarks with the finite values of the current masses. The 



calculated neutron-proton mass difference appeared to be strongly density dependent. The result 
overestimated the value of NSA. The application of the bag models were considered by Hatsuda et al. 



134 1 . The chiral bag model provided the proper sign of the effect, but underestimated its magnitude. 



5.3.2 QCD sum rules view 

The QCD sum rules look to be a reasonable tool for the calculation of neutron-proton binding energy 
difference for the nucleons, placed into the isotope-symmetric nuclear matter. Denoting 

5x = x n — x p (191) 

for the strong interaction contribution to the neutron-proton difference of any parameter x, we 
present 

5e = 5U + 5T (192) 

with U(T) — the potential (kinetic) energy of the nucleon. To the lowest order in the powers of 
density one finds 

ST = - % 5m , (193) 
m 2 

while the value 5U = 5m m can be obtained from the sum rules. 

The attempts to apply QCD sum rules for solving the NSA problem were made in several papers 



134| | [ p~3T| . We shall follow the papers of Drukarev and Ryskin ||137|| which present the direct 



extension of the approach, discussed above. It is based on Eqs. (|182|) - (|184j) with the terms £™ being 



calculated with the account of the finite values of the current quark masses. Besides the quark mass 
difference the CSB effects manifest themselves through isospin breaking condensates 

_ (0\dd-uu\0) (M\dd-uu\M) 
70 ~ (0|H0> ; 7m ~ (M\uu\M) ' ( } 

The characteristics 70 and 7 m are not independent degrees of freedom. They turn to zero, at 
\i = 0, being certain (unknown) functions of //. The dependence 70 (a0, 7m(/-0 can be obtained in 
framework of the specific models. Anyway, due to the small values of m u ^ we expect I70I, \^ m \ *C 1. 
Following the strategy and keeping only the leading terms, which are linear in /1, we shall obtain the 
energy shift in the form 

5e = din + a 2 7o + a 3 7 m (195) 
with a,i being the functions of the density p (the contribution 5T is included into ai). 

To calculate the mass-dependent terms in lhs of Eqs. ( P-82j) -( I51 ) one should include the quark 



masses to the in- medium quark propagator — Eq. ( |149| ) . The first term in rhs of Eq. (|149|) , which is 
just the free quark propagator should be modified into {ix a p) / {2n 2 x A ) — m q / (2tt 2 x 2 ). This provides 
the contribution to the scalar structure of the correlator in the lowest order of OPE. Account of the 
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finite quark masses in the second term of Eq. (|149 ) manifest themselves in the next to leading orders 
of OPE. Say, during evaluation of the second term in Subsection 5.2, we used the relations expressed 
by Eq . (|53"D , which were obtained for the massless quarks. Now the first of them takes the form 

1 Til ■ — 

(#> = i (<M-^(N\m\N) (196) 

for the flavour % = u,d. This leads to the contribution to the vector structure of the correlator pro- 
portional to the scalar condensate. Also, the second moment of the scalar distribution is proportional 
to the vector condensate, contributing to the scalar structure of the correlator. 

The leading contribution caused by the scalar structure of the correlator 

(5U) 1 = 0.18 ^ ^ GeV , (197) 
m p 

while the CSB term, originated by the vector structure is 

(SU) 2 = - 0.031 * K{p) ~ K{po) GeV (198) 
m p 

with po — 0.17 Fm" 3 being the value of saturation density. The term A£™ of Eq.( |186| ) provides the 
contribution containing the CSB condensate 7 m — Eq. ( |194|) 

(5U) 3 = 32 ( 7r >(p) - «(0)) + ( 7 m - 7o)«(0)) / GeV 2 . (199) 



For the complete calculation one needs the isospin breaking shifts of vacuum parameters 5X 2 and 
5W 2 while the empirical value of 5m can be used. Thus, the analysis of CSB effects in vacuum 
should be carried out in framework of the method as well. This was done by Adami et al. [ |135|| . The 
values of 5\ 2 , 5W 2 and of the vacuum isospin breaking value 70 were obtained through the quark 
mass difference and the empirical values of the shifts of the baryon masses. This prompts another 
form of Eq.(TO 

5s = b l p + b 2 - fm (200) 

with b\ = a\ + a 2 ( 7 o//i), b 2 = a 3 . The expressions for the contributions to 5e, caused by the 
shifts of the vacuum values 5m, 5\ 2 and 5W 2 are rather complicated. At p = p the corresponding 
contribution is 

(5U) 4 = -0.4 MeV . (201) 
For the sum J2i(5U)i we find, after adding the contribution 5T 

61(A)) = -0-73 , b 2 {p ) = -1.0 GeV . (202) 



The numerical results can be obtained if the value of 7m is calculated. This can be done in 
framework of certain models. However, even now we can make some conclusions. If we expect the 
increasing restoration of the isospin symmetry with growing density, it is reasonable to assume that 
\lm\ < I To I - Also, all the model calculations provide 70 < 0. Thus we expect 70 < 7m < 0. If 7 m = 
we find 5e = —2.4 MeV, eliminating the vacuum value 5m = 2.3 MeV. Hence, the isospin invariance 
appears to be restored for both the condensates and nucleon masses. 
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The present analysis enables also to clarify the role of the CSB effects in the scalar channel. 
Indeed, neglecting these effects, i.e. putting {8U)2 = (SU)s = we obtain 5e > 0. This contradicts 
both experimental values and general theoretical expectations. Thus we came to the importance of 
CSB effects in the scalar channel. 

Adami and Brown ||I35|| used NJL model, combined with BR1 scaling for calculation of parameter 
7 m . They found 7 m /7o = ( K (p)/ K (0)) 1//3 - Substituting this value into Eq. fl200| ) we find 

6s = (-0.9 ± 0.6)MeV (203) 

with the errors caused mostly by uncertainties of the value of 70. A more rapid decrease of the ratio 
7m/7o would lead to larger values \5e\ with 5e < 0. Putting 7 m = 70 provides Se = — 0.3MeV. 

Of course, Eq. ( |203| ) is obtained for infinite nuclear matter and it is not clear, if it can be ex- 
trapolated for the case A = 40. We can state that at least qualitative explanation of the NSA is 
achieved. 



5.3.3 New knowledge 

As we have stated earlier, the QCD sum rules can be viewed as a connection between exchange 
of uncorrelated qq pairs between our probe nucleon and the matter and the exchange by strongly 
correlated pairs with the same quantum numbers (the mesons). In the conventional QHD picture 
this means that in the Dirac equation for the nucleon in the nuclear matter 

(q-V)tp = (m + $)^ (204) 

the vector interaction V corresponds to exchange by the vector mesons with the matter while the 
scalar interaction $ is caused by the scalar mesons exchange. In the mean field approximation the 
vector interaction V is proportional to density p, while the scalar interaction is proportional to the 
"scalar density" 

f d 3 p m* 

Wp^iM- (205) 

which is a more complicated function of density p — see Eqs. (|76|) , (|78|) . Thus V = V(p), while 
$ = $(p s )- We have seen that QCD sum rules provide similar picture in the lowest orders of OPE: 
vector and scalar parts of the correlator G m depend on vector and scalar condensates correspondingly: 
G™ p = G™ p {v(p))] G™ = G™(/t(p)). However as we have seen in Subsection 5.2.6, we find a somewhat 
more complicated dependence in the higher order OPE terms, say, G™ = G™(K,(p),v(p)), depending 
on both scalar and vector condensates. This means that the corresponding scalar interaction $ = 
$(p s ,p), requiring analysis beyond the mean field approximation. 

As one can see from Eqs. (|197|) and (|198|) in the case of CSB interactions such complications 
emerge in the sum rules approach in the leading orders of OPE. 

Thus the QCD sum rules motivated CSB nuclear forces V and $ in Eq.( |204| ) are expected to 
contain dependence on both "vector" and "scalar" densities, i.e. V = V(p,p s ) and $ = $(p, p s ). As 
we said above, such potentials can emerge due to complicated structure of nucleon-meson vertices. 
This can provide the guide-lines for building up the CSB nucleon-nucleon potentials. 
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Another new point is the importance of the CSB in the scalar channel. Neglecting the scalar 
channel CSB interactions we obtain the wrong sign of the effects, i.e. Se > 0. This contradicts the 
earlier belief that the vector channel u — p mixing is the main mechanism of the effect ||132| . Our 



result is supported by the analysis of Hatsuda et al. ||139|| who found that the uo — p mixing changes 
sign for the off-shell mesons. This can also help in constructing the CSB nuclear forces. 



5.4 EMC effect 



The experiments carried out by EMC collaboration ||140|| demonstrated that deep inelastic scattering 



function F^(xb) of nucleus with atomic number A ( Xb stands for Bjorken variable) differs from the 
sum of those of free nucleons. Most of the data were obtained for iron (Fe). The structure function 
was compared to that of deuteron, which imitates the system of free nucleons. The deviation of the 
ratio 

^ (lB ) = (206) 

from unity is caused by deviation of a nucleus from the system of free nucleons. The ratio R{xb) 
appeared to be the function of xb indeed. Exceeding unity at xb < 0.2 it drops at larger xb reaching 
the minimum value R Fb (xb) ~ 0.85 at xb ~ 0.7. This behaviour of the ratio was called the EMC 
effect. 

There are several mechanisms which may cause the deviation of the ratio R(xb) from unity. These 
are the contribution of quark-ant iquark pairs, hidden in pions, originated by the nucleon-nucleon 
interactions, possible formation of multiquark clusters inside the nucleus, etc. Here we shall try to 
find how the difference of the quark distributions inside the in-medium and free nucleons changes 
the ratio R(xb)- 

The QCD sum rules method was applied to investigation of the proton deep inelastic structure 
functions in vacuum in several papers. The second moments of the structure functions were obtained 



by Kolesnichenko ||141|| and by Belyaev and Block | 142 |. The structure function ^(xg) at moderate 



values of xb was calculated by Belyaev and Ioffe [122]. Here we shall rely on the approach, developed 



by Braun et al. ||143|| which can be generalized for the case of finite densities in a natural way. On 
the other hand, such generalization is the extension of the approach discussed in this section. 

To obtain the structure function of the proton, the authors of [[L43|| considered the correlation 



function G, describing the system with the quantum numbers of proton, interacting twice with a 
strongly virtual hard photon 

G(q,k) = i 2 J d 4 xd 4 y e t{qx)+i{ky) (0\T[ V (x)f](0)]H(y, A)\0) . (207) 

Here q and q + k are the momenta carried by the correlator in initial and final states, k = k\ — k 2 is the 
momentum transferred by the photon scattering. The incoming (outgoing) photon carries momentum 
kiik?), interacting with the correlator in the point y — A/2 (y + A/2). The quark-photon interaction 
is presented by the function H(y,A). In the next step the double dispersion relation in variables 
q\ = q 2 and q\ = (q + k) 2 is considered. The crucial point is the operator expansion in terms of the 
nonlocal operators depending on the light-like (A 2 = 0) vector A [ 143 1 . After the Borel transform 



in both q\ and q 2 is carried out and the equal Borel masses M\ = M 2 are considered, the Fourier 
transform in A provides the momentum distribution of the quarks. 
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Figure 14: Second order interaction of the hard photon (dashed line) with the correlator. The other 
notations are the same as in the previous pictures. 
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Figure 15: The in-medium changes of the d quark distribution (dashed curve) and of the u quark 
distribution (dash-dotted curve) of the fraction x of the momentum of the target nucleon. The solid 
curve presents the function R — 1 with the ratio R, defined by Eq.(206). 
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This approach was applied by Drukarev and Ryskin ||144| for calculation of the quark distributions 
in the proton, placed into the nuclear matter. The two types of contributions to the correlator should 
be considered — Fig.l4a,b. In the diagram of Fig. 14a the photon interacts with the quark of the 
free loop. In the diagram of Fig. 14b it interacts with the quark exchanging with the matter. The 
modification of the distributions of the quarks was expressed through the vector condensate, which 
vanishes in vacuum and through the in-medium shifts of the other condensates and of the nucleon 
parameters m, A 2 and W 2 . The result appeared to be less sensitive to the value of the four-quark 
condensate than the characteristics of the nucleon considered in Subsection 5.2. 

Note that the results are true for the moderate values of x only and cannot be extended to the 
region x C 1. This is because the OPE diverges at small x 
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Omitting the details of calculation, provided in ||144|| we present the results in Fig. 15. One can 
see that the distributions of u and d quarks in fraction of the target momentum x are modified in a 
different way and there is no common scale. The fraction of the momentum carried by the u quarks 
(x u ) decreases by about 4%. The ratio R, determined by Eq.( |206| ) has a typical EMC shape. 



The technique used in [144] can be expanded for the calculation of the quark distributions at 
1 < x < 2. Thus the approach enables to describe the cumulative aspects of the problem as well. 



5.5 The difficulties 

In spite of the relative success, described above, the approach faces a number of difficulties. Some of 
them take place in vacuum as well. The other ones emerge in the case of the finite density. 

The first problem is the convergence of OPE in the lhs of the sum rules. Fortunately, the 
condensates which contribute to the lowest order OPE terms can be either calculated in a model- 
independent way, or expressed through the observables. This is true for both vacuum and nuclear 
matter — at least, for the values of density which are close to the saturation values po- However 
the situation is not so simple for the higher order OPE contributions. The four-quark condensate 
is the well known headache of all the QCD sum rules practitioners. The problem becomes more 
complicated at finite density, since the conventional form of presentation of this condensate contains 
the strongly cancelling contributions. 

In order to include the higher OPE terms one needs the additional model assumptions. The same 
is true for the attempts to go beyond the gas approximation at finite densities. Recently Kisslinger 



145|| suggested a hybrid of QCD sum rules and of the cloudy bag model. 

Note, however that QCD sum rules is not a universal tool and there are the cases when OPE does 
not converge. We mentioned earlier that this takes place for the nucleon structure functions at small 
x - [ 122 1 . Some time ago Eletsky and Ioffe ||147|| adduced the case when the short distance physics 
plays important role, making the OPE convergence assumption less convincing. Recently Dutt- 
Mazumder et al. [|146|| faced the situation when the ratio of two successive terms of q~ 2 expansion is 



not quenched. 

There are also problems with the rhs of the sum rules. The "pole+continuum" model is a very 
simple ansatz, and it may appear to be oversimplified even in vacuum. The spectrum of the nucleon 
correlator in-medium is much more complicated than in vacuum. The problem is to separate the 
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singularities of the correlator, connected with the nucleon from those of the medium itself. As we 
have seen, the "pole+continuum" model can be justified to the same extent as in vacuum until we 
do not include the three-nucleon interactions. We do not have a simple and convincing model of the 
spectrum which would include such interactions. 



Anyway, the success of vacuum sum rules [117] and reasonable results for the nucleons at finite 
densities described in this Section prompt that the further development of the approach is worth 
while. 



6 A possible scenario 

The shape of the density dependence of the quark scalar condensate in the baryon matter n(p) 
appears to be very important for hadronic physics. It is the characteristics of the matter as a whole, 
describing the degree of restoration of the chiral symmetry with growing density. On the other hand, 
the dependence n(p) is believed to determine the change of the nucleon effective mass m*(p). The 
shape of the dependence m*(n(p)) differs in the different models. 

The lowest order density dependence term in the expansion of the function n(p) is model- 
independent. However for the rigorous calculation of the higher order terms one needs to know the 
density dependence of the hadron parameters m*(p), m* A (p), f*(p), etc. In Sec. 4 we presented the 
results of the calculations of the condensate n(p) under certain model assumptions. A more detailed 
analysis requires the investigation of the dependence of these parameters on QCD condensates. 

Such dependence can be obtained by using the approach, based on the in-medium QCD sum rules. 
In Sec. 5 we show how in- medium QCD sum rules for the nucleons work. Even in a somewhat skeptical 



review of Leinweber ||148|| , where the present state of art of applications of the QCD sum rules is 
criticised, the method is referred to as "the best fundamentally based approach for investigations of 
hadrons in nuclear matter" . Of course, to proceed further one must try to overcome the difficulties, 
discussed in Subsec.5.5. 

The lowest order condensates can be either calculated, or connected directly to the observables. 
This is true for both vacuum and nuclear mater. However, neither in vacuum nor in medium the 
higher order condensates can be obtained without applications of certain models. Thus in further 
steps we shall need a composition of QCD sum rules with model assumptions. 

We have seen that the density dependence of the delta isobar effective mass is important for the 
calculation of the nonlinear contribution to the scalar condensate n(p). The shape of this dependence 
is still obscure. Thus the extension of the QCD sum rules method for the description of A-isobars in- 
medium dynamics is needed. Such work is going on — see, e.g., the paper of Johnson and Kisslinger 
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The fundamental in-medium Goldberger-Treiman and Gell-Mann-Oakes-Renner relations are 
expected to be the other milestones of the approach. The agreement with the results with those of 
conventional nuclear physics at p ~ p would be the test of the approach. 

Further development of the approach would require inclusion of the vector mesons. The vector 
meson physics at finite densities is widely studied nowadays. Say, various aspects of QCD sum rules 
application where considered in recent papers ||146| , |150| , |151| , |152|| while the earlier works are cited 
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m reviews 



F21, PS 



We expect the investigation in framework of this scenario to clarify the features of baryon pa- 
rameters and of the condensates in nuclear matter. 
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